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C»H. Abstract. We present an outlook for possible discovery of supersymmetry with broken R-parity at Run II 

^ ■ of the Tevatron. We first present a review of the literature and an update of the experimental bounds. In 

1-^ ' turn we then discuss the following processes: 

1. Resonant slepton production followed by |?p decay, (a) via LQD'^ and (b) via LLE'^ . 
• T-H , 2. How to distinguish resonant slepton production from Z' or W' production. 

ji^ ■ 3. Resonant slepton production followed by the decay to neutralino LSP, which decays via LQW^. 

^ ' 4. Resonant stop production followed by the decay to a chargino, which cascades to the neutralino LSP. 

C^ I 5. Gluino pair production followed by the cascade decay to charm squarks which decay directly via LiQ2D^. 

6. Squark pair production followed by the cascade decay to the neutralino LSP which decays via L1Q2D1. 

7. MSSM pair production followed by the cascade decay to the LSP which decays (a) via LLE'^, (b) via 
LQD", and (c) via WWD", respectively. 

8. Top quark and top squark decays in spontaneous ftp. 



I INTRODUCTION 

H. Dreiner 

The chiral superfields of the supersymmetric Standard Model are 

i. (1,2,-i), E^ {1,1,1), i/i(l,2,-i), if2(l,2,i), 

Q, (3,2,^), U^{3,l,-h, D'r(3,l,\). (1) 

6 6 6 



In parentheses we have given the quantum numbers under the Standard Model gauge group, Gsm = SU{'i) x 
SU{2) X C/(l). I = 1,2,3 is a generation index. The most general superpotential for this field content is^ 

W^WMSSM + Wfi^, (2) 

Wfi^ = ^X^.keabLtL'^.Ei + K^^^ea,L1Q]Dl + X'l^^e^y^Ur Df Dl^ + n^e^b^Hl (3) 

Here, Wmssm is the superpotential of the MSSM containing the terms which give mass to the SM fermions. 
We do not further discuss it here. i,j,k = 1,2,3 are generation indices, a,b ~ 1,2 are SU{2) isospin indices, 
and x,y, z = 1, 2, 3 are SU{3) colour indices. A, A', A" are dimensionless Yukawa couplings. Ki are mass terms 
mixing the Icptonic and Higgs doublets. The terms in Wfi violate the discrete and multiplicative symmetetry 
R-parity 

where B is baryon-, L is lepton- number and S denotes the spin. It is the purpose of this chapter to investigate 
the phenomenological consequences of the superpotential terms Wf^ at the Tevatron. For a recent review on 
Rp violation (_^p) see [1]. A very nice detailed overview of |?p collider phenomenology is given in [2]. 

At the level of the superpotential, the terms LiH2 can be rotated away by a field redefinition [3]. When 
including the soft-breaking terms this is no longer true in general. However, we shall split our investigations 
into explicit Rp violation, where we neglect all contributions from the Lii?2 mixing, and spontaneous ^p. 

II EXPLICIT R-PARITY VIOLATION 

A Existing Bounds 

B. Allanach, A. Dedes, H. Dreiner 

In Table 1, we show a recent compilation of indirect bounds on the ^p Yukawa couplings [4,5,1,6]. We present 
the bounds at the 2a level and give two significant figures where possible. The bounds depend on the mass 
of super partner and we give the functional dependence. For more detail on how the bounds are obtained see 
[4,5]. We ignore strict bounds from cosmology [7] because they are model dependent [8]. A discussion of the 
boimds on A" j, is given in Sect. VIII. 

B Outline of Signals and Summary of the Literature 

The two main differences between the standard MSSM phenomenology and the phenomenology of explicit 
^p are 

1. Single production of supersymmetric particles is possible. For example at a hadron collider we can 
have resonant slcpton (charged and neutral) production via the operators LiQjDk and resonant squark 
production via the operators UiDjDk- The lowest order Feynman diagrams are shown in Figs. la. This 
lowers the kinematic threshold for the discovery of supersymmetry. 

2. The lightest supersymmetric particle (LSP) is not stable, and can possibly decay in the detector. If the 
LSP is the lightest neutralino it can decay, for example, via the operator LiQjDk, as shown in Fig. 2. 
However, since the LSP is not stable, the cosmological argument requiring it to be electrically and colour 
neutral no longer applies [20] . To date there have been no systematic studies of non- neutralino LSPs in 
the context of ^p. If the LSP is charged and has an appreciable lifetime (or is stable), it can be detected 
in a search for charged massive particles (CHAMPs). 

A first systematic study of R-parity violating signals at hadron collider was presented in [21,22]. There are in 
principle two kinds of signals given by the two kinds of supersymmetric production. 



Note that '^ here refers to charge conjugation and is not a 5(7(2) index. 



TABLE 1. Latest 2g limits on the magnitudes of weak scale trilinear R-parity 
violating couplings from indirect decays and perturbativity. We have made use 
of the data in [17]. The dependence on the relevant superparticle mass is shown 
explicitly. When the perturbativity bounds are more stringent than the empirical 
bounds for masses mj -=1 TeV, then we display them in parentheses. Where a 
bound without parentheses has no explicit mass dependence shown, the mass de- 
pendence was too complicated to detail here and a degenerate sparticle spectrum 
of 100 GeV is assumed. 

ijk X^jk{Mzr X,fc(Mz)'^ K'AMzY 

111 - 5.2 X IQ-* X™ 

112 - 0.021 x,,;^ io-^x(^^f/2 

113 - 0.021 x^jji^ 10-* 

121 0M9xj^^^ 0.043 x^^jj^^ Iq-^^ ^ (^^)V2 

122 0.049 X ,„,/",, 0.043 x 



100 GcV 'J.u-t^ ^ lOj^GcV 

123 0.049 x,^;^ OMSXj,^ (1.23) 

131 0.062 x^j^;;^ 0.019 x^jjo^ 10"* 

132 0.062 X j^;^^ 0.28 x j^^^^ %^ (1.04) (1.23) 

133 O.OOeOv'j^/if/lOO GeV 1.4 x IQ-^Jm-jWO GeV 
211 0.049 X ,„„"", 0.059 x 



100 GcV "-"^-^ ^ 100 GeV 

212 0.049 x,^^ii^ 0.059 x,^^^ (1.23) 

213 0.049 X ,5^;^ 0.059 X ,^5^±g^ (1.23) 

221 - 0.18 X j^^ (1.12) (1.23) 

222 - 0.21x^^^^(1.12) 

223 - 0.21 X j^^ (1.12) (1.23) 

231 omOXj^^ 0.18 X ^j3j^±|^ (1.12) (1.23) 

232 0.070 X wgIv 0. 56 (1.04) (1.23) 

233 0.070 X loo^lv O.lSv'rn.s/lOO GeV 

311 0.062 x^^|i^ o.llx^^(L12) 

312 0M2Xj^^^^_ 0.11x^^^^(1.12) 0.50(1.00) 

313 O.OOOOymf/lOO GoV 0.11 x ww (1-12) 0.50 (1.00) 

321 0.070 X j;^^ 0.52 x ^^^^ (1.12) 0.50 (1.00) 

322 0.070 X ,j3j^4^ 0.52x^^^(1.12) 

323 0.070 X j^^^ 0.52 x ^^j^ (1.12) 0.50 (1.00) 

331 - 0.45 (1.04) 0.50 (1.00) 

332 - 0.45 (1.04) 0.50 (1.00) 

333 ; 0.45 (1.04) - 

^ Updated bounds from Rof. [5,1]. Bounds on A121, A122, A123 have been obtained 
from charged current universality [9]. Bounds on A131, A132, A231, A232 and A233 
have been derived from [9] measurements of Rt = r(r -^ evh')/V{T -^ ^vu) 
and Rrp, = V{t —> iivv)/V{ii —> euv) [17]. The bound on A133 [10] has been 
obtained from the experimental limit on the electron neutrino mass [17]. 

''Bounds on Aii2, Aii3, Ai2i, Ai22, and Ai23 have been obtained from charged 
current universality [9]. The bound on \'\ii has been derived from neutrino- 
less double beta decay [11,12] where ™ = (me'/lOO GeV)^ x (m^o/100 GeV)'''^ 
and on Ai3i from atomic parity violation [9,15]. This latter bound is at the 
3(T level, since the data disagree with the standard model at the 2.5cr level 
[15]. The bound on Ai32 comes from the forward-backward asymmetry in e+e~ 
collisions [9]. Bounds on A'133, A233 have been obtained from bounds on the 
neutrino masses [10] and on A211, A212, A213 from i?^ = r(7r -^ ev)/T{Ti -^ ^v) 
[9,5]. Bounds on A221, A231 come from v^ deep inelastic scattering [9,5] and 
on A222, A223 from the D-meson decays [5], D ^> Klv. The bounds without 
parentheses on A232, A331, A332, A333 have been derived from Ri = T{Z -^ 
had)/r(Z -^ If) for m, = 100 GeV [13] and on Agn, A312, A313 from Rr-n = 
r(r -^ 'Ki/T)/T{n —> nvii) [9,5]. The bounds on the couplings A321, A322 and A323 
have been derived from Ds decays [5], i.e., Rd^ — r{Ds — > ri^T-)/r(Ds — + /ii^^). 
There are also bounds on A3J3 from Rt [16] but these are weak at 2a level and 
thus not displayed. 

'^ The indirect bounds on X'/j/. existing in the literature are on A112 from double 
nucleon decay [14] (A is a hadronic scale and it can be varied from 0.003 to 
1 GeV and {j^f^^ from 2 x lO" to 10^ for m,- = 100 GeV) and on A113 
from neutron oscillations [18,14] for m,-=100 GeV. For m5=200 (600) GeV 
the bound on A113 is 0.002 (0.1). Bounds on X'^j^. have been derived from 
Ri = r{Z -^ had)/r(Z -^ II) at la for m = 100 GeV [19] and, for heavy 
SQuark masses, is not more strineent than the Derturbativitv bound, which is 
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FIGURE 1. Production of x°^+ 






FIGURE 2. 

as well. 



LQD decay of x'^. The neutralino is a Majorana fermion and decays to the charge conjugate final state 



1. Single Sparticle Production, e.g. u + d ^ f+. 

2. Pair production of Sparticles via MSSM couplings, e.g. u + u ^ q + q. 

In both cases the production is followed either by direct R-parity violating decays of the supersynimetric 
particles or by cascade decays to the LSP, which in turn decays via ^p. This latter decay can be prompt 
or delayed. If it is delayed it can lead to a detached vertex or decay outside the detector. To date only 
prompt decays have been considered. Within this report we also discuss delayed decays (see VIII and [23]) 
leading to detached vertices. The LSP decays outside the detector for small R-parity violating couplings, 
A, A', A" ~ 10~^, [22,24]. The resonant production is then also suppressed and we thus retrieve the MSSM 
phenomenology, where sparticles are pair-produced and the LSP escapes detection. This is discussed in the 
contribution by the SUGRA working group. 

In the following we summarize the possible signatures and indicate what work has already been performed. 
There is at present only one published experimental search for ^p at hadron colliders [25] . A summary as well as 
the projection of this analysis for Run II is presented in Sect. VI. There is a preliminary DO analysis publically 
available [26]. We present the extension of this analysis to Run II in Sect. VII. There are published results on 
searches for pair production of leptoquarks, which can be interpreted as searches for ^p (see Sect.IIB2). In 
the following, most references are thus to theoretical work on specific signals. 

1 Single Sparticle Production 

At a hadron collider sneutrinos and charged leptons can be produced on resonance via the operator LiQjD'^. 
They can then decay again via the same operator, another LQD'^ operator, the operator LiLjE'^ or via the 
MSSM gaugino interactions. Thus we obtain the possible reactions 



djdk 



dkdj, LiQjDk, {b) 
ly^xL MSSM, (c) 



Xr, 



MSSM, (d) 



ijak 



i^Vj, ULjEk, (e) 

Ujdk, LiQjDk, (/) 

i+xl, MSSM, (g) 

ly^xt., MSSM. 



(5) 



(h) 



Here we have indicated the decay products of the slepton, as well as the coupling which leads to this decay. 
The decays to gauginos proceed via MSSM gauge couplings. Resonant slepton production was first considered 
in [21], where the production cross sections for the Tevatron were determined and the decays (5)a,b,e,f were 
compared to the background. The specific case of scalar tau production in (5)e was considered in [27] and 
compared to more recent Tevatron data on Drell-Yan di-lepton production. In Sect. Ill we investigate the 



reactions (5)a,b,e,f in detail and compare them to exisiting Tevatron data. Furthermore the reactions (5)a,e 
can look identical in their final state to new gauge boson, W' , Z' , production and we also investigate how this 
can be disentangled in Sect. III. For both cases we use present data to project the search reach for Run II. In 
Sect. IV, we study for the first time the case of like-sign di-lepton production via the reactions (5)g. The focus 
is on a single dominant ftp operator, so that the neutralino decays as x° ^^ ^^ + 2 jets, as shown in Fig. 2. The 
analysis for the reaction (5)d is completely analogous. 

The reactions for resonant squark production via the operators U'^D'^D'^ are 



an' -^ n"^ i ^^ ' u u ly , \^u) , . 





' (^,i^)»'Z, 


L,QD^, 


(a) 


n"^ 


qq', 


WD^D", 


ib) 


q ^ 


q{x"m,9), 


MSSM, 


(c) 




^ Qxti^ 


MSSM. 


(d) 



The two-jet production via (6)b was considered in [21]. In special cases the reactions (5)f and (6)b can lead 
to single top quark production. This has been studied in [28]. In Sect.V we study the reaction (6)d in detail, 
where the produced squark is a stop. On resonance, the reaction (6)c is kinematically blocked because of the 
final state top quark. The chargino of reaction (6)d cascade decays to the lightest neutralino in a semi-leptonic 
decay and the charged lepton is used for detection. The dominant decay of the neutrlino is also suppressed 
because of the top quark in the final state and thus decays outside the detector. 

2 Sparticle Pair Production 

For ^p couplings A, A', A" ^ 10^^ resonance production is no longer viable. We then expect pair production 
of supersymmetric particles to dominate. The main reactions are 

pp-" {99, m, m, x°x°, xtx] , x°x^, i'^r, w, ^+i>} (7) 

The supersymmetric particles can undergo cascade decays to the LSP as in the MSSM. The LSP, which we 
here assume to be the lightest neutralino then decays via ^p 






The resulting reactions were first considered in [21] and classified in [22]. Gluino and squark pair production 
were first analysed in detail and compared to CDF data in [29]. In each case, the MSSM cascade decays were 
explicitly neglected, giving branching ratio 1 for the decay to the neutralino LSP. In [30] the production of the 
electroweak gauginos was studied focusing on the reactions 

"^^ X?X? - (X?^±^)X?, (c) ^'^ 

I X^xr ^ (X?^+^)(X?^'-^), id) 

followed by the decay of the LSP via the operator LiLjE^ as in (8). In this case the discovery reach is 
significantly above that of the MSSM. 

In Ref. [31] all the production mechanisms given in Eq.(7) were studied simultaneously in the framework 
of minimal SUGRA to see what combined contribution there would be to several possible discovery signals: 
(i) quartic-leptons, (ii) tri-leptons, (iii) opposite sign di-leptons, (iv) same-sign di-leptons, and (v) missing 
transverse energy. Each signature was compared for the two scenarios of an LSP decaying via LLE'^ and 
WD'^D'^, respectively. In Sect. VI we study the case of pair production followed by cascade decay to the LSP 
including a detector simulation. The LSP decays via the operator LLE'^. The results are compared to the 
present Tevatron data and projected to a discovery reach for Run II. In Sect. VII we investigate in detail the 
potential at Run II for the di-lepton signature in the case of a dominant LQD'^ operator. In [32] the previous 
analysis of gluino and squark pair production of [29] was extended to include cascade decays and investigate the 



^Ut^,^ 


ULjEl, 


(a) 


£f + 2 jets. 


UQjDl 


(b) 


I' + 2 jets. 


UQjDl, 


(c) 


3 jets. 


U^D'^PI. 


id) 



search reach in the gluino and squark mass at the Tevatron. In Sect. VIII, we investigate the possibihty of the 
LSP decaying via VD'^D'^ to three jets, where the chargino and ncutraUno mass are nearly degenerate, thus 
making the cascade leptons of [31] too soft for detection. This is the worst case scenario for supersymmetry 
searches. 

In February, 1997 both HERA experiments announced an anomaly in their high Q^ deep inelastic scattering 
data [33]. This has not been confirmed in the later data but it has also not been excluded yet. The initial 
anomaly could naturally be explained in terms of resonant production of squarks via ^p [34] or as resonant 
production of leptoquarks [35] . It lead to an increased interest in the searches for squarks and leptoquarks at 
the Tevatron. In leptoquark pair production at the Tevatron, the leptoquarks directly decay to a lepton and 
a quark. The experimental signature is equivalent to pair produced squarks decaying directly via the LQD^ 
operators. This is thus a distinct signature from the cascade decay to the LSP. We can interpret the leptoquark 
searches [36,37] in terms of ^p with large AJ j, Yukawa couplings. We would expect additional decays for the 
squarks to gauginos. So the leptoquark mass bound corresponds to an excluded region in the mass branching 
ratio, BR{q -^ £ + q'), plane. The present Tevatron data [36-38] excludes the solution of the HERA anomaly 
where BR{q ^ £ + q') = 1 and severely constrains the case where BR{q -^ £ + q') < 1 [39]. In Sect. VI we 
present a study of gluino pair production followed by the cascade decay to charm squarks (as motivated by 
the HERA anomaly) which in turn decay directly via L1Q2D1, i.e. not via the LSP. We also investigate the 
pair production of squarks followed by the cascade decay to the LSP which decays via L1Q2DI. This focuses 
particularly on the ^p interpretation of the HERA data which is not covered by the leptoquark searches. 

PART 1: Resonant SParticle Production 

III R-PARITY VIOLATING DECAY OF THE SLEPTON AND 

IDENTIFICATION 

J. Hewett, T. Rizzo^ 

In this section, we concentrate on the two sets of trilinear i-violating terms in Wj^, , and the case of single 

charged or neutral slepton production via qq^ ' annihilation at the Tevatron through the A' couplings. If this 
slepton decays to opposite sign leptons (through the A couplings) then an event excess, clustered in mass, will 
be observed in the Drell-Yan channel similar to that expected for a new neutral or charged gauge boson, Z' 
or W'. In addition, both £ and D resonances may decay hadronically via the same vertices that produced 
them, leading to potentially observable peaks in the dijet invariant mass distribution. Thus resonant slepton 
production, first discussed in Ref. [21,40], offers a unique way to explore the ftp model parameter space at 
hadron colliders. It is important that ^p also allows for other SUSY particles, such as t and/or 6, to be 
exchanged in the non-resonant t, u— channels and also contribute to Drell-Yan events. However, it can be easily 
shown that their influence on cross sections and various distributions will be quite small if the low energy 
constraints on the Yukawa couplings are satisfied. [35] The questions addressed in this analysis are: (i) what 
are the mass and coupling reaches for slepton resonance searches at the Tevatron in the Drell-Yan and dijet 
channels and (ii) how can slepton resonances, once discovered, be distinguished from Z' , W production. 

3 Drell- Yan Channel for Sleptons 

From the bounds in Section II A, we see that the coupling responsible for first generation slepton production, 
A']^]^]^, is much too restricted to generate a sizable cross section. However the bounds on the couplings for second 
and third generation slepton production and decay are all of similar numerical values and allow for a reasonably 
sized signal. In addition, the background free e^/i^ signature is also possible. 

In the case of Drell-Yan production the search reach analysis is straightforward being nearly identical to 
that used for new gauge boson production, apart from acceptance issues due to the differences between spin-0 
and spin-1 resonances. Since sleptons are expected to be narrow, the narrow width approximation is adequate 
and the analysis presented in Ref. [41] can be directly followed. In addition to the slepton mass the only other 
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FIGURE 3. Discovery regions (lying below the curves) in the mass-coupling plane for Jfip resonances in the neutral 
(left) and charged (right) Drell-Yan channels at the Run II Tevatron. From top to bottom the curves correspond to 
integrated luminosities of 30, 10, 5 and 2 fb~^. The estimated reach for Run I is given by the lowest curve. The 
parameter X is defined in the text. 

parameter in the calculation is the product of the appropriate Yukawa couplings. A', from the initial state du or 
dd coupling vertex, and the slepton's leptonic branching fraction, Bi. Denoting this product as X = (A')^i?£, 
we obtain the 95% C.L. search reach as a function of X in both the charged and neutral channels. These 
results are displayed in Fig. 3 for various values of the integrated luminosity. Not only is it important to notice 
the very large mass reach of these colliders for sizeable values of X ~ 10"'^, but we should also observe the 
small X reach, X ~ lO^'^'"'^^-' and below, for relatively small slepton masses. These results show the rather 
wide opportunity available to discover slepton resonances over extended ranges of masses and couplings at the 
Tevatron. Note that for fixed values of X the search reach is greater in the charged current channel due to the 
higher parton luminosities in that case. 





m, (GeV) m, (GeV) 

FIGURE 4. Cross sections for narrow dijet resonances (solid) at the 2 TeV Tevatron arising from v (left) or £ (right) 
production in comparison to the anticipated search reaches of CDF (dots) and DO (dashes). The upper(lower) curve for 
each experiment assumes an integrated luminosity of 2(30) fb^^. The three solid curves from top to bottom correspond 
to slepton resonance predictions for Y—Q.l, 0.01 and 0.001, respectively, where Y is defined in the text. 



4 Di-Jet Decays of the Slepton 



Since dd and/or ud annihilation are responsible for single slepton production, slepton resonance can also 
decay into these same fermion pairs. Hence ^ or ;> will have dijet decays which may appear as observable 
peaks above the conventional QCD backgrounds. Searches for such narrow dijet resonances have already 
been performed [42] at the Tevatron by both CDF and DO during Run L Using their results and scaling by 
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FIGURE 5. Comparison of the Monte Carlo generated normalized angular distribution for the leptons in Z' decay 
(data points) with that for a O (solid curve) assuming a (left)400 event sample or a (right)2000 event sample; the 
displayed errors are purely statistical. The Z' is assumed to have no forward-backward asymmetry due to its fermionic 
couplings. A 1 + z^ angular distribution is represented by the dotted curve. 

appropriate factors of beam energy and integrated luminosities we may estimate the probable search reaches 
for CDF and DO from Run II. Our estimates conform to the expectations given in Ref. [43]. The cross sections 
themselves are calculable in the narrow width approximation in terms of the product Y = (A')^i?2j, where A' 
is the Yukawa coupling and B2J is the dijet branching fraction. The results are presented in Fig. 4. Here we 

see that for values of F ~ 0.001 — 0.01 the Tevatron will have a substantial mass reach for slepton induced dijet 
mass bumps during Run II. Note that as in the case of Droll- Yan production, larger cross sections for fixed 
Y occur in the charge current channel than in the neutral current case due to the larger parton luminosities. 
Unfortunately, if such a bump is observed it will not be straightforward to identify it as a slepton resonance. 



5 Distinguishing Sleptons from W', Z' 

Our next issue concerns the ability to determine that an observed resonance (or Jacobian peak) in the Drell- 
Yan channel is due to slepton production, instead of a new gauge boson, Z' , W . One clear signature for v 
production would be the observation of e/x final states. These are not expected to occur in Z' models and 
would point to ^p. In addition, if the Rp conserving decay modes of the slepton dominate there will be no 
identification problem. If the ]^p modes dominate one should first look for universality violations, e.g., if the 
resonance decays to only one of e+e^ or fi^fi^ or if these two rates are substantially different. Most models 
with new gauge bosons do not have substantially different couplings to the first two fermion generations. 

In the neutral current channel, i.e. , i> or Z' production, the forward-backward asymmetry Afb in the leptonic 

decay distributions provides a good diagnostic tool. It is well-known that most models of Z' bosons have parity 

violating fermionic couplings which would lead to Aps 7^ 0. However, i>, being spin-0, would always produce a 

null asymmetry. Aps is more easily measured and requires less statistical power than does the reconstruction 

of the complete angular distribution. This is important since, whereas only 10 or so background free events 

would constitute a discovery, many more events, ~ 100 — 200 are required to determine the asymmetry. This 

implies that the reach for performing this test is somewhat, if not substantially, less than the discovery reach. 

For example, the Tevatron may discover a i) with a mass of 700 GeV for a certain value of X but only for 

masses below 500 GeV would there be enough statistics to extract Afb for this same X value. 

A more complex and interesting situation arises when the Z' naturally has Apg := as in, e.g. , some Eq 

models [44]; in this case the on-resonance asymmetry data alone is insufficient. If Apg could be measured 
throughout the resonance region, it would be possible to deduce through detailed line-shape studies whether 
or not the new contribution interferes with the SM amplitude (something that does not occur in the case of 
spin-0 i> production). However, finite dilepton mass resolution, especially for the /i^/i~ final state, in addition 
to the increased required statistics may disrupt this program. 

With a plethora of statistics the complete angular distribution can be obtained as shown in Fig. 5. Here 
we compare Monte Carlo generated data for a Z' with a zero forward-backward asymmetry with both the flat 
distribution expected in the spin-0 case, and the ^ 1 + z^ distribution (z ~ cos 9) of spin-1 exchange. We 





FIGURE 6. The lepton charge asymmetry in the charged current DreU-Yan production channel at the 2 TeV Tevatron 
for the SM (sohd curves) and with 250(700) GeV £ exchange (the dash-dotted and dashed curves) assuming A, A' = 0.15. 
From top to bottom in the center of the figure, the SM curves correspond to Mt bins of 50-100, 100-200, 200-400 and 
> 400 GeV, respectively. 



ignore complications due to possible acceptance losses arising from rapidity cuts in the forward and backward 
directions. Such a distribution has been measured by CDF both on and above the SM Z resonance [45]. 
This analysis indicates that of order ~ 1000 events are required to make a clean measurement, which is a 
sample approximately 100 times larger than that required for discovery. Although such measurements would 
be conclusive as to the identity of the spin of the resonance, the required statistics results in a significant loss 
in the mass range over which it can be performed. In the example discussed above where the slepton search 
reach was 700 GeV we find that the angular distribution could only be determined for masses below ^ 400 
GeV assuming the same X value. 

Angular distributions cannot be used to separate between possible £ and W' production in the charged 
current channel due to the missing energy in the event. However, there are two other useful observables in 
this case. First, is the transverse mass (Mt) distribution associated with the new Jacobian peak region where 
interference with the SM amplitude would occur for W' , but not slepton, production. This test is unfortunately 
far more difficult to perform than in the neutral current channel again due to the missing energy as well as 
mass smearing. A second possibility is to examine the leptonic charge asymmetry, ^(77), for the electrons or 
niuons in the final state as a function of their rapidity, where 



A{v) 



dN+/dr]^dN^/dr] 
dN+/dr] + dN^/dr] 



(10) 



where N± is the number of positively/negatively charged electrons of a given rapidity, t]. In the SM, the charge 
asymmetry is sensitive to the ratio of u-quark to d-quark parton densities and the V ~ A nature of the W 
production and decay. [46] Since the coupling structure of the SM W has been well-measured, any deviations in 
this asymmetry within the Mt bin surrounding the W mass have been attributed to modifications in the parton 
distributions. Here, we are more interested in events with larger values of Mt- Note that A{—ri) — —A{ri) if 
CP is conserved (which we assume) so that we will only need to deal with 77 > in the following discussion. 

Figure 6 displays the lepton charge asymmetry within four Mt bins, corresponding to 50 < Mt < 100 GeV, 
100 < Mt < 200 GeV, 200 < Mt < 400 GeV and 400 < Mt < 1800 GeV, for the SM and with a 250(700) 
GeV £ with, for purposes of demonstration. A, A' = 0.15. In particular, the lepton charge asymmetry can be 
significantly altered for larger values of Mt in the bins associated with the new Jacobian peak. Note, however, 
that there is essentially no deviation in the asymmetry in the transverse mass bin associated with the SM W 
peak, 50 < Mt < 100 GeV, so that this Mt region can still be used for determination of the parton densities. 
The figure also shows that the presence of the £ tends to drive the asymmetry to smaller absolute values as 
perhaps might be expected due to the presence of a spin-0 resonance. 

Figure 7 presents the corresponding modifications in the leptonic charge asymmetry due to an 800 GeV W 
with either purely lcft-handed(LH) or purely right-handcd(RH) fcrmionic couplings. Note that the W with 
purely RH couplings, unlike the LH W' , does not interfere with the SM amplitude, similar to the case of £ 
production. We see that the deviation in the asymmetry due to either type of W is very different than that 
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FIGURE 7. Same as the previous figure but now for the case of a 800 GeV W' with purely left-handed(left panel) or 
purely right-handed(right panel) couplings. 
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FIGURE 8. Direct comparison of the charge asymmetry induced by a 800 GeV ^(solid) and a left- or right-handed 
W^'(dot and dash) of the same mass, taking 600 < Mr < 900 GeV. The assumed values of the Yukawa couplings are as 
in the earlier figure. 

for a slepton. In the W' case, the magnitude of the asymmetry is substantially increased for both left- and 
right-handed couplings and these two cases are themselves potentially distinguishable by using the data in the 
Mt bin below, but not containing, the Jacobian peak. 

The Mt bins we have taken in this analysis are rather broad. We might expect that if the width of the 
Mt bin around the W' or slepton Jacobian peak is compressed then the purity of the resonant contribution 
will be increased resulting in a better separation of the two cases (at the price of reduced statistics) . These 
expectations are realized in Fig. 8 which shows a more direct comparison of the lepton charge asymmetries for 
a i and W of the same mass (800 GeV) taking the width of the Mt bin surrounding the Jacobian peak to be 
only 300 GeV. However, the left- and right-handed W possibilities are no longer separable since this requires 
interference with the SM. 

Deviations in the leptonic charge asymmetry can also be used to probe indirectly for the exchange of i 
through ftp couplings. To demonstrate this, we fix the £ width to mass ratio to be T/m — 0.004 and subdivide 
each of the four Mt bins discussed above into rapidity intervals of A77 ==0.1. For a fixed slepton mass and 
integrated luminosity we generate Monte Carlo data for various values of the Yukawa couplings and perform a 
X^ analysis to obtain the sensitivity to the product AA'. The 95% C.L. search reach from this analysis is shown 
in Fig. 9 where it is clear that the reach obtained in this manner is rather modest. 

6 Summary 



The above analysis shows that resonant s-channel slepton production with subsequent decay into purely 
leptonic or dijet final states via^^p couplings is observable in hadronic collisions over a wide range of parameters. 
If this signature is observed, we have demonstrated that the leptonic angular distributions and the lepton charge 



0.100 
0.050 



0.010 
0.005 



5 10 15 20 25 30 

L (tb-') 
FIGURE 9. Search reach for i exchange as a function of the integrated luminosity assuming V jvn = 0.004 for masses 
of 1500, 1000, 750 and 250 GeV (from top to bottom). 



TABLE 2. 


Summary of the Background Simulation 


Process 


Cross section 
before Cuts /nb 


No. of Events Expected No. of Events after cuts, 
simulated for 2fb^^ luminosity 


bb mixing 

tt 

single top 


(9.3 ±2.3) X 10"* 
6.81 ±0.31 
1.55 ±0.12 


8.3x10® 0.12 ±0.12 
2.0x10^ 0.02 ±0.02 
3.5x10'* 0.00 ±0.03 


Total 


0.14 ±0.13 



asymmetry can be successfully used to distinguish slcpton resonances from those associated with new gauge 
bosons. This process provides a clean and powerful probe of _^p supersymmetric parameter space. 

IV R-PARITY CONSERVING DECAY OF THE SLEPTON 

R. Dreiner, M. Seymour, P. Richardson 

In this section we discuss the resonant charged slcpton production via the operator LiQjDf. followed by the 
decay to the neutralino decay as shown in (5)g. 



J-jCLk 



11 



it + 



An 



The neutralino in turn then decays via the same _^p-operator LjQjDk, 



x° 



{£■ +uj+dk; t'i±dj±dk}. 



(11) 



(12) 



Since the neutralino is a Majorana fermion it can also decay to the charge conjugate final states, with equal 
probability. In spirit, this is similar to the HERA process considered in [47,48]. 

The tree-level Feynman diagrams for the slepton production and the neutralino decay are shown in Fig. 1 
and Fig. 2, respectively. Due to the Majorana nature of the neutralino, we have a signature of two like-sign 
charged leptons. In the following we shall consider only electrons or muons, i.e. we focus on the operators 
LcQjDk and L^QjDk- We expect these leptons to have high transverse momentum, px, and be well isolated 
whereas the leptons from the Standard Model backgrounds have lower px and arc also poorly isolated. We 
therefore hope that this signature can be seen above the background if we apply isolation and px cuts. 



A Backgrounds 

In the following we combine the backgrounds for both electrons and muons. The main backgrounds to this 
like-sign dilepton signature are as follows 



Fig 3a. Both leptons with P^ above 15 GeV/c 
800 



Fig 3b. Both leptons with P^ above 20 GeV/c 
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FIGURE 10. Effect of the Parton-Level Cuts. 



1. bb production followed by the production of at least one B^ ^ meson, which undergoes mixing. If the two 
b-quarks in the event decay semi-lcptonically this gives two like-sign charged leptons. 

2. tt production followed by t -^ W+b -^ c^iJ^h, and t -^ W^b -^ qqb -^ qqW+c — > qqe+i^cC. 

3. Single top production (s and t channel) followed by semi-leptonic decays of the top and the B-meson 
produced after hadronization. 

4. Non-physics backgrounds from fake leptons and charge misidentification. There are also backgrounds 
due to the production of weak boson pairs, i.e. WZ and ZZ, where at least one of the charged leptons is 
not detected [49]. These require a full simulation including the detector. We do not consider them here. 

We use HERWIG 6.0, [50-52], to simulate these background processes. The program includes the compu- 
tation of the supersymmetric spectrum and the MSSM decay branching ratios from the ISASUSY program 
[53] . Due to the high cross section for the production of bb it was necessary to impose a parton- level cut of 
20GeV on the px of the b and b to enable us to simulate a sufficient number of events. In Fig. 10 we show 
the distribution of events (using the full Monte Carlo simulation) as a function of the (parton-level) px of the 
bottom quark for two different values of the lepton px cut. We did not simulate any events for which the px of 
the bottom quark was below 20 GeV since the cross section is too large. If we extrapolate using Figs. 10a, b to 
lower b-quark px we can see that for a lepton px cut of 20 GeV, Fig. 10b, our approximation should be good, 
i.e. we expect the area under the curve for px(b) < 20 GeV to be neghgible. For px(^) > 15 GcV, Fig. 10a, we 
would still expect a significant number of events at 15 GeV < px(b) < 20 GcV. Besides the parton-level cut, 
we forced the B-mesons to decay semi-leptonically. This means we neglect the production of leptons from the 
decay of charmed mesons which should also be a good approximation as we expect the leptons produced from 
these decays to be poorly isolated. 

Figure 11 displays the effect of the lepton isolation cut on the bb background for two different values. The 
effect of the isolation cut on the tt and single top backgrounds is shown in Fig. 12. As can be seen in Figs. 11, 12, 
by imposing an isolation cut of 5 GeV and a cut on the px of the leptons of 20 GcV the background can be 
almost eliminated. Table 2 shows the backgrounds with a px cut on the leptons of 20 GeV and an isolation 
cut of 5 GeV. We have used the leading-order cross section for the bb and single top backgrounds and the 
next-to-leading order cross section, with next-to-leading-log resummation, from [54] for the tt cross section. 
In both cases the error on the cross section is the effect of varying the scale between half and twice the hard 
scale, and the error on the number of events is the error in the cross section and the statistical error from the 
simulation added in quadrature. Realistically we cannot reduce these statistical errors due to the large number 
of events we would need to simulate. We have implemented the full hadronization using HERWIG 6.0. 

With these cuts and using Poisson statistics, a 5a fluctuation of the total background corresponds to 4 events 
with an integrated luminosity of 2fb^^. Hence we consider 4 signal events to be sufficient for a discovery of 
the new ^p signal process. 
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Fig 4a. Effect of an Isolation Cut of 3 GeV 
on the bottom background 
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Fig 4b. Effect of an Isolation Cut of 5 GeV 
on the bottom background 
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FIGURE 11. Effect of the lepton isolation cut on the bb background. The upper curve is the full bb background and 
the lower curve is obtained after imposing the isolation cut. 



B Signal 



To simulate the signal and the effect of the cuts, we modified HERWIG 6.0, [50-52], to include the production 
process, the MSSM decay of the slepton, and the ^p decay of the neutralino. The decay rate of the neutralino 
and its branching ratios were calculated in the code and a matrix element for the neutralino decay [55,56] was 
implemented in the Monte Carlo simulation. 

We use the program to estimate the acceptance of the signal process, i.e. the fraction of the like-sign dilepton 
events which pass the cuts multiplied by the branching ratio to give a like-sign dilepton event. Fig. 13 shows 
the acceptance for two different SUGRA points, with an isolation cut on the leptons of 5 GeV and a cut 
Pt(^) > 20 GeV. As can be seen in Fig. 13b, the acceptance drops in two regions. For lower values of Mq, the 
slepton is not much heavier than the neutralino. The charged lepton from the decay of the slepton is then quite 
soft and gets rejected by the px cut. For large values of Mq the slepton is much heavier than the neutralino. 
The neutralino then gets a significant boost from the slepton decay. The neutralino decay products are folded 
forward in the direction of this boost causing the event to be rejected by the lepton isolation cut. 

To estimate the acceptance properly we need to run a scan of the SUGRA parameter space using the Monte 
Carlo event generator. This still remains to be done. To give some idea of what range of couplings and masses 
we may be able to probe instead we assume an acceptance of 10% using the same cuts as before. We can then 
estimate the range of couplings which may be accessible. 

As can be seen in Figs. 14, 15 the production cross section for A211 — 10^^ is sufficient to produce a signal 
which is more than 5cr above the background for large regions of the SUGRA parameter space. In some regions 
where the neutralinos become Higgsino-likc the cross section drops. The cross section also drops as we approach 
the region where the neutralino is heavier than the slepton and the resonance becomes inaccessible. 

We focused on the coupling A211 because the experimental bound on X'm from neutrinoless double beta 
decay is very strict, [57,1], as discussed in Sect. II A. The bound on X'ln weakens as the squark mass squared 
and for squark masses above about 300 GeV (which we expect in the SUGRA scenario for the heavier slepton 
masses) X'm w 10^^ is experimentally allowed and our analysis thus applies for this case as well. A211 ~ 10~^ 
is well within the present experimental bounds [1]. 

In these figures we see that we are sensitive to slepton masses up to 300 GeV for couplings of 10~^. The 
production cross section scales with the square of the coupling. For slepton masses around 100 GeV, just above 
the LEP limits, we can thus probe couplings down to about 2 x 10^'^. 



Fig 5a. Effect of an Isolation Cut of 5 GeV 
on the top background 
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Fig 5b. Effect of an Isolation Cut of 5 GeV 
on the single top background 
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FIGURE 12. Effect of the Isolation Cut on the tt and single top Backgrounds 
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FIGURE 13. Acceptance from Monte Carlo Simulation for two different SUGRA points. 



C Conclusion 



We have performed an analysis of the physics baekground for hke-sign dilepton production at Run II and 
find that with an integrated luminosity of 2 fb^^, a cut on the transverse momentum of the leptons of 20 GeV 
and an isolation cut of 5 GeV the background is 0.14 ± 0.13 events. This means that 4 signal events would 
correspond to a 5ct discovery, although in a full experimental analysis the non-physics backgrounds must also 
be considered. 

Using a full Monte Carlo simulation of the signal including a calculation of the neutralino decay rate, its 
partial widths and a matrix element in the simulation of the decay we found that the acceptance for the signal 
varies but for a reasonable range of parameter space is 10% or greater. 

When we then look at the cross section for the production of x^^"*" we find that we can probe ^p couplings 
of 2 X 10^^ for slepton mass of 100 GeV and up to slepton masses of 300 GeV for ^p couplings of 10^^, and 
higher masses if the coupling is larger. 



Acknowledgments: The work of P.R. is supported by PPARC, UK. 






c 
o 

u 

0) 
CO 



o 

u 
o 



M2=100GeV/c'= 



tan|8 = 2 



DiscoVery limit ' ' 

- Slepton Mass 100 GeV/c^ 

\ - - - Slepton Mass 300 GeV/c'^ I 


- 


1 



/x/GeV/c" 



^ 

d 



C 
o 

o 
m 



o 
o 



M2 = r)00GcV/c^ tari.S = 2 



- — Discovery limit 

Slepton Mass 100 GeV/c^ 

. _ Slopton Mass 300 GeV/c' 




200 400 

/i./GeV/c^ 



FIGURE 14. Cross section for tan/3 — 2 and A211 — 10^'^. The long-dashed horizontal line corresponds to the 
discovery limit of 4 signal events for 2fb~^. The signal cross section after cuts is given by the solid and dot-dashed 
curves for two slepton masses. 






a 
o 

u 

0) 

m 



O 

o 



M2=100GeV/c'' tan,e=l( 



- — — Discovery limit 

Slepton Mass 100 GeV/c^ 

- - - Slepton Mass 300 GeV/c^ 






o 

o 



o 



M2 = 500GeV/c^ tanjS=10 



— — — Discovery limit 

Slepton Mass 100 GeV/c^ 

Slepton Mass 300 GeV/c* 



Ai/GeV/c^ 

FIGURE 15. Cross section for tan/3 = 10 and Aan = 10"^ 



/Li/GeY/c' 



V SINGLE-TOP-SQUARK PRODUCTION VIA UDD COUPLINGS 

E.L. Berger, B. W. Harris, and Z. Sullivan 



In supersymmctric extensions of the standard model, the bounds on possible i?-parity-violating couplings 
are relatively restrictive for the first two generations of quarks and leptons, but much less so for states of 
the third generation (c/. Sect. II A. If Rp is conserved superpartners must be produced in pairs. In many 
models, the squarks and gluinos are relatively heavy, and therefore their pair production incurs a large phase 
space suppression at the energy of the Tevatron. Here, the s-channel production of a single squark through an 
-^p-mechanism [21] is considered. Our motivation is that the greater phase space offsets the reduced coupling 
strength in the production. We focus on the relatively light top squark ti [58] and its subsequent Rp- conserving 
decays [59]. The Rp violation penalty is paid only once, in the initial production, and is offset by the greater 
phase space relative to pair production. 

We choose one clean i?p-conserving decay mode, ti -^ bxi ^ with Xi ^ l'^ + u + x\- Here, I is an electron 
or muon, and the Xi ^^d Xi ^I'e the chargino and lowest-mass neutralino states of the MSSM. For top-squark 
masses in the range of 180--325 GeV, we simulate both the signal and standard model background processes 
and thereby show that the top squark can be discovered, or the current bound on the size of the _Rp-violating 
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FIGURE 16. Cross section for 7?p-violating production of a single top-squark at Run II of the Tevatron {X^ij — 0.1) 
compared with the cross section for i?p-conserving production of pairs. Both cross sections are computed at leading 
order. 

coupling Ag'j ■ can be reduced by up to one order of magnitude with existing data and by two orders of magnitude 
at the forthcoming run II of the Fermilab Tevatron [59] . 

In general it is possible to have _Rp-violating contributions to the MSSM superpotential of the baryon- or 
lepton-number violating type. However, limits on the proton decay rate severely restrict their simultaneous 
presence. We therefore assume the existence of a baryon-number-violating coupling only of the form [60] 



Wi?p = ^.kU^DWl 



(13) 



Here, U^ and D'^ are right-handed-quark singlet chiral superfields, i,j, k are generation indices, and c denotes 
charge conjugation. For production of a right-handed top squark via an s-channel diagram d^d'^ — > u}^ the 
relevant couplings are A3j^2i -^3131 ^^'^ -^323 • The current best limits on these couplings come from experiments 
at the CERN Large Electron Positron Collider (LEP). Measurement of Ri, the partial decay width to hadrons 
over the partial decay width to leptons of the Z boson, limits the values of AJj'j. to be less than approximately 
0.5 (c/Sect. HA). 




FIGURE 17. Branching ratio for the top squark to decay into two jets via the i?-parity-violating coupling A" as a 
function of the coupling. The values of the SUGRA parameters are given in the text. 



The colour- and spin-averaged partonic cross section for qq' -^ ti is 



TABLE 3. Cuts used to simulate the ac- 
ceptance of the detector at the Tevatron 
run II, and run I (in parentheses if differ- 
ent). The lepton veto (Eti < 45 GeV) is 
optimized for small top-squark mass. 



I'7fl<2(l) 
|?7i|<2.5(l.l) 
|»7>|<2.5 
\ARij\ >0.7 
fij, > 20 GcV 



ETb > 40 GcV 

Eti > 15 GeV (20 GeV) 

Etj > 20 GeV 

\AR,i\ >0.7 
Eti < 45 GeV 



27r 
"3 



^~l-^3y 



, sm 



2 

ti 



-6(1- ml /s) 



(14) 



where \/s is the partonic center of mass energy, and 9^ relates the right- and left-handed top-squark interaction 
states to the mass eigenstates. The hadronic cross section depends on the following combinations of incident 
parton distribution functions (PDF's): d (8) s, d (E) b, and s (E) b, where d, s, and b denote the PDF's of the 
down, strange, and bottom quarks, respectively. We assume Agj^j ~ -^313 ~ -^323 — ^" ■ ^^^^ numerical results 
represent the sum of ti and ti production. 

The mass dependences of the cross sections for single and pair production of top squarks at the Fermilab 
Tevatron {y/s = 2 TeV) differ significantly, as shown in Fig. 16. The curves are based on CTEQ4L leading 
order parton distribution functions [61] and A" ~ 0.1, an order of magnitude below the current limit. In fact, 
A" can be smaller by an additional order of magnitude and the single-top-squark cross section will still exceed 
that for pair production if mj > 100 GeV. 
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FIGURE 18. The reconstructed- mass M distribution for single-top-squark production plus background (S + B) at the 
Tevatron (y's = 2 TeV) for a top-squark mass mj_^ = 242 GeV. Also shown is a decomposition of the total background 
(B) into it's major components, all W jet modes (Wj) and single-top-quark production (t). The coupling A" = 0.03 
produces the minimum signal for a 5a significance at this mass. 



In the MSSM, the up-type squark u*^ can decay into charginos and neutralinos via the two-body processes 



1,2,3,4). Three body modes are possible, including 



h 



-^ d^ +XJ (.?■ = 1,2) and u'' 
W^ + fe + Xij which is similar to decay into the top quark (followed by top decay to VF+ -(- 6) but softer; 



«'= + X^ (.7 



and ii -^ c -|- x? via a flavor-changing loop process. The squark can also decay into a gluino plus a quark if 
allowed kinematically. 



With baryon- number- violating terms present, the top squark can also decay into quark pairs via the A" 
couplings. The branching fraction into two jets is shown in Fig. 17 for the specific SUGRA point discussed 
below. The branching ratio is quite insensitive to the SUGRA parameters over a wide range. If A" is large, the 
decay to quark jets dominates. However, as shown below, the i?p-conserving decay still produces a measurable 
and useful cross section. 

For the remainder of this study, we focus on the two-body decay mode [22] ti — > b + Xi , with x+ -^ l + u + Xi- 
Here, I denotes an electron or muon, which usually comes from a W. We expect Xi to undergo an _Rp-violating 
decay outside of the detector, so it is treated as stable for our purposes. 

A minimal supergravity model [62] is adopted to obtain the relevant masses and decay branching fractions. 
We begin with common scalar and fermion masses of toq — 100 GeV and TOi /2 = 150 GeV, respectively, at the 
GUT scale. We choose a trilinear coupling Aq — —300 GeV and the ratio of the Higgs vacuum expectation 
values tan/3 = 4. The absolute value of the Higgs mass parameter /i is fixed by electroweak symmetry breaking 
and is assumed positive. Superpartncr masses and decay widths are calculated with IS A JET [53] . At the weak 
scale, m^ = 183 GeV, m^o == 55 GcV, and rn~± — 103 GeV. In order to isolate the effects of the ^p sector, 
we vary mo and keep the other supersymmetric parameters fixed. Since the gaugino masses depend primarily 
on the choice of rni/2, variation of toq allows us to vary m^ without any appreciable change in the masses of 
the decay products. 
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FIGURE 19. Statistical significance of the single-top-squark signal {S/vE) in run II of the Tevatron (y^ = 2 TeV, 
2 fb^^) versus A" for a variety of top-squark masses. 



For the four-body decay mode of interest, the signal consists of a tagged 6-quark jet, a Icpton, and missing 
transverse energy associated with the unobserved neutrino and Xa. The dominant backgrounds, in order of 
importance, arise from production and decay of the standard model processes Wc, in which a lU is produced 
in association with a charm quark c that is mistaken for a fe; Wj, in which a lU is produced in association 
with a hadronic jet that mimics a 5; Whb; Wcc; and single-top-quark production via Wg fusion. For these 
background processes, we work with tree-level matrix elements obtained from MADGRAPH [63] convolved 
with leading-order CTEQ4L [61] parton distribution functions, at a hard scattering scale fi^ = s. In an 
experimental analysis, the Wj background will be normalized by the data. To simulate the resolution of the 
hadron calorimeter, we smear the jet energies with a Gaussian whose width is AEj/Ej = 0.80/-^^© 0.05 
(added in quadrature). 

We simulate the acceptance of the detector by using the selections listed in Table 3. The assumed coverage 
in rapidity for taggable &-quark jets and leptons is smaller for run I than for run II. However, the signal and 
background are similar in shape in these variables, and thus S/ B is not sensitive to this cut. The lepton must 
be isolated from any jets, as defined by a cone of radius AR, or it is considered missed. Similarly, an isolation 
cut is used for the 6-quark jet in order to help identify it. We assume a 6-tagging efficiency of 60% (50% 
for run I) with a mistag rate of 15% for charm quarks and 0.5% for light quarks [43]. Missing momentum is 
identified with a massless particle whose three-momentum balances that of the 6-quark and the Icpton. 

As expected from the primary decay, ti ^ b + xt , the distribution in the transverse energy Et of the b 



quark is peaked sharply near the maximum value allowed kinematically. The spectrum of the background b 
quark is soft, and thus we impose a hard cut {ETb > 40 GeV) on the minimum Et of the fe-quark jet. For 
larger values of mj this cut would be optimized by raising the minimum Exb- The fo-jet becomes too soft to 
be detected if imq < m^-^ + E'™*. This contributes to a lower limit on mj below which our proposed search 
mode is not useful. 

The background from single-top- gwarfc production will produce a peak in any mass reconstruction. We utilize 
the fact that single top-quarks are often produced with extra hard jets, and we impose a "jet veto" . Wc require 
that there be no hard jets in the hadron calorimeter {Etj > 20 GeV, |?7| < 2.5), beyond the one that is 6-tagged. 
After the jet veto, the remaining background is due almost entirely to misidentified charm and light-quark jets 
from Wc and Wj production. The transverse energy of the lepton tends to be relatively soft for the signal 
at low m^ , whereas it peaks around 40 GeV when it comes from the W in the background. A cut to remove 
hard leptons, with Eti > 45 GeV, reduces the background by a factor of 2 with little effect on the signal at 
low masses. The final significance for the signal at run II is barely changed by this "lepton veto", but it is 
especially helpful for the run I data. 
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FIGURE 20. Lower limits on discovery (S/VB = 5), evidence {S/vB = 3), and 95% confidence-level exclusion 
{S/yB — 1.96) for A" versus top-squark mass in Run I of the Tevatron {^/s — 1.8 ToV, 110 pb^^), and in Run II 
(V^ = 2TeV, 2fb-i). 



Shown in Fig. 18 is an example of the signal and background. For this case, mj ~ 242 GeV and A" ~ 0.03. 
The mass variable is defined as M^ = (Pf, -I- P; -|- PxY where the Pb and Pi are the four- momenta of the h and 
lepton. The four-momentum Px is defined such that its three-momentum balances that of the h and lepton, 
and P|- = 0. The clearly discernible signal in Fig. 18 would constitute a discovery at the level of 5ct with an 
integrated luminosity of 2 fb"-'^ at ^/s = 2 TeV. When mq is reduced to 183 GeV, the signal and background 
spectra peak at about the same location, and sensitivity to the signal begins to be lost. 

Examination of the structure of the cross section involving Rp- conserving decay modes reveals that as A" 
grows, the decrease in branching fraction is compensated by the increase in cross section. As depicted in 
Fig. 19, 
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where f{Rp) is a constant times the branching fraction into i?p-conserving modes. As A" — > oo, the cross 
section goes to a constant; whereas, when A" -^ the cross section decreases as |A"p. The relationship 
S/^/B (X \X"\'^/\/B, valid for small A", implies a lower limit on the values of A" that can be probed. On the 
other hand, this relationship highlights an insensitivity to variations in the estimate of the background. 

In Fig. 20, we show the reach in A" for 180 < mj < 325 GeV. With an integrated luminosity of 2 fb^^ at 
vo = 2 TeV, discovery at the level of 5a is possible provided that A" > 0.02-0.05. Otherwise, a 95% confidence 



level exclusion can be set for A" > 0.01-0.03. For the lower integrated luminosity and energy of the existing 
Run I data, values of A" > 0.03-0.2 can be excluded at the 95% confidence level if m^ =180-285 GeV. 

We conclude that, as long as the i?p-conserving decay ti -^ bxi -^ it^Xi is allowed, it should be possible to 
discover the top squark or to reduce the lower limit on A" by two orders of magnitude at Run II of the Fermilab 
Tevatron, for 180 < m^^ < 325 GeV and Ag^^ > 0.02-0.05. Existing data from Rim I of the Tcvatron should 
allow a reduction of the limit on A" by an order of magnitude. With such a reduction, one can establish that 
i?p-violating decay is unlikely and rule out most of the possible influence of the top squark on single-top-quark 
production and decay. 

This work was supported by the U. S. Department of Energy, High Energy Physics Division, under Contract 
No. W-31-109-Eng-38. 

PART 2: SParticle Pair Production 

VI GLUINO AND SQUARK PAIR PRODUCTION WITH LEPTONIC 

DECAYS 

M. Chertok 

CDF has performed searches for squark and gluino production followed by -^p-decays using data taken with 
the Run I detector [64] . Upgrades of both the detector and Tevatron are currently underway [65] . These will 
provide substantial enhancements for these and other searches in Run II, scheduled to begin in 2000. During 
this run, CDF II will collect roughly 2 fb~ of data at y/s — 2 TeV, corresponding to twenty times the present 
statistics. The 10% increase in energy corresponds to a 40% increase in the ti yield, and similarly will aid new 
phenomena searches. 

A ^p Decays via LQD'^: Di-Lepton Events 

Events with a positron and a jet at high Q^ values, detected at the HERA experiments [33], have sparked 
interest in ^p since [34] such events can be explained by the production and decay of a charm squark (c^): 
e^ + d ^f cl -^ e^ -\- d, where Rp is violated at both vertices. For this scenario, c^ with mass M{cl) — 
200 GeV/c^ is the preferred squark flavor, because its associated coupling A']^2i [66,67] is less constrained by 
experiment than the others. 

We test two |?p processes that involve the same A']^2i coupling [25]: 

(16) 
(17) 

The first we denote the "charm squark analysis" and the second the "neutralino analysis". Here, the |?p decays 
are indicated by "=>." For process (16) we assume M{g) > M{cl) = 200 GeV/c^. The masses of the other 
squarks are given in a MSSM scenario in Rcf. [66]. For process (17), we consider qq production (5 degenerate 

squark flavors) and iiii production separately. We also make the mass assumptions: M(xi ) > M{q) > 
M(xi), M{xf) « 2M(xi), and M{xf) > M{ii) - M{b), where the first relation suppresses q -^ x^, the 
second relation arises from gaugino mass unification, and the third ensures that Br{ii -^ ex?) = 100% for 
M(ti) < M{t). Given the Majorana nature of the gluino and neutralino, reactions (16) and (17) yield like-sign 
(LS) and opposite-sign (OS) dielectrons with equal probability. Since LS dilepton events have the benefit of 
small SM backgrounds, we search for events with LS electrons plus two or more jets. 

This analysis requires at least two electrons with Et > 15 GeV in the central electromagnetic calorimeter. 
The r]-(j) separation ARee = ■\/(A(/i)2 -f- (A?])^ between two electrons must be greater than 0.4. Each electron 
must pass an isolation cut which requires the total calorimeter Et in an ?]-((> cone of radius AR ~ 0.4 around 
the electron, excluding the electron Et, to be less than 4 GeV. Jets are identified in the calorimeter using cone 
size AR = 0.7 for clustering and in the range |?7j| < 2.4. There must be at least two jets with Et > 15 GeV, 
ARjj > 0.7, and AR^j > 0.7. We further require no significant fk in the event: prl \/Y1,^t < 5 GeV^/^. No 
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FIGURE 21. Exclusion regions in the M{q)-M{g) plane (left) and upper limits on cross section times branching ratio 
for fiti production (right). 



LS ee events survive our selection. The remaining 166 OS ee events are consistent with expected contributions 
from SM processes, notably Drell-Yan production of dielectrons. 

Event acceptances are calculated using Monte Carlo samples generated with IS A JET, CTEQ3L parton 
distribution functions, and passed through the CDF detector simulation program. SM backgrounds for this 
search, tt and bb/cc production, are small. We find the total background in 107 pb~ is consistent with zero 
events. We set limits on the cross section times branching ratio for the two processes under study. For the 
charm squark analysis, the event acceptance is a very weak function of M{g) in the range of 16.0% to 16.6%; 
we include a 10% systematic uncertainty (dominated by the uncertainty on the integrated luminosity), and 
exclude a ■ Br > 0.18 pb independent of M{g). In Figure 21 we plot the results from the charm squark analysis 
in the gluino-squark mass plane. Contours are shown for two values of the branching ratio Br{cL —* ed), where 
we have compared our results to the NLO gg production cross section multiplied by the branching ratio to LS 
ee from Ref. [67]. 

For the neutralino analysis, the acceptance is determined for each squark and neutralino mass pair and 
ranges from 3.7% to 15.2%. We obtain a ■ Br limits which range as a function of the squark mass from 0.81 
pb to 0.26 pb for a light neutralino, and from 0.35 pb to 0.20 pb for a heavy neutralino. Figure 21 shows 
the results for the neutralino analysis for the case of iiii production. Plotted are our 95% C.L. upper limits 
along with the cross section times branching ratio versus M{ti) from the NLO prediction. The branching ratio 
Br{ti -^ cxi) is taken to be 1.0 We also assume Br{x'i -^ qq'e) — Br{xi — ^ qq't^) — 1/2, which need not be 
the case [47]. Since each neutralino decays to e"*" or e~ with equal probability, the branching ratio to LS ee 
is 1/8. This analysis excludes M{ii) below 120 (135) GeV/c^ for a light (heavy) neutralino. The results for 
the neutralino analysis for the case of five degenerate qq production are displayed in Figure 22. Again, plotted 
is our cross section times branching ratio limit for two neutralino masses, along with the NLO prediction 
which includes a gluino mass dependent i-channel contribution to the cross section. Thus, we set gluino and 
neutralino mass-dependent lower limits on the degenerate squark mass in the range from 200 to 260 GeV/c^. 

We include a Run II projection for this search in Figure 22 which folds in the increased cross section at 
^/s — 2 TeV along with the ratio of the integrated luminosities. Assuming event selection and background 
levels similar to those in Run I, this search will be sensitive to squark masses up to 380 GeV/c^. The addition 
of the muon channel will further enhance the reach, as will the improved lepton identification in Run II. 



B 



Decays via LLE^: Multilepton Events 



If we now assume one of the LLE'^ operators is non-zero, the LSP decays as in Eq.(8)a. For the case of 
supersymmetric pair production as in Eq.(7), followed by the cascade decay to the LSP, we then expect events 
with extra leptons. CDF has an analysis underway that examines A121 ^ [68]. As in [31], all production 
processes (except scalar top) are considered simultaneously. The coupling A121 allows two classes of processes: 
1) direct decays of sleptons and sneutrinos to leptons and neutrinos, and 2), decays of charginos and ncutralinos 
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FIGURE 22. Upper limits on cross section times branching ratio for five-flavor degenerate qq production (left) and 
projection for Run II (right). Also shown are the NLO cross sections multiplied by a branching ratio to LS ee of 1/8 
with the assumption of Br{q -^ gXi) — l-O- 

TABLE 4. Multilepton ftp search projections for Run II. The 
efficiencies are based on the Run I detector simulation. 
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to leptons and neutrinos as in (8)a. It is further assumed that A121 is large enough so that the ^p decays are 
prompt, i.e. within 0(lmm) of the interaction point, but also that A121 <C 1. This tends to suppress both 1) and 
2), except for decays of the LSP which have no competing MSSM decays. Thus, Xi ~^ Mci'e, fieve, eev^ eev^ 
each with equal probability. Such events therefore contain four (or more) leptons and f}^. 

Simulation of the signal is performed using ISAJET within a SUGRA framework. The mass parameters Mq 
and Ml 12 are input, while the others are fixed: Aq — 0, tan/3 = 2, sgn/i = —1. The event selection requires 
four leptons, £ = (e,/x), where one lepton satisfies Erie) > 12GeV, \i]e\ < 1-1 ot prifj) > 12 GeV/c, \r]fj,\ < 0.6, 
while the others must satisfy Erie) > 5GeV, \r]e\ < 1.1 or pt(m) > 5 GeV/c, |ryp| < 1.0. Furthermore, the 
leptons must be separated by AR > 0.4. Neither lepton isolation nor ^ requirements are imposed. SM 
processes that can yield four lepton events include tt, bb/cc, and Z'^Z'^ production and decay, although the 
rates are very low. Another background under study is that of processes containing three real leptons plus an 
additional fake lepton. 

Projections for this search in Run II are presented in Table 4, which folds in the approximate factor of twenty 
in statistics and the increase in ^/s. The background increase to 0(10) events in Run II could be controlled 
by additional cuts on ^, jets, or lepton isolation. Improvements to the detector, for example the upgraded 
plug EM calorimeter, will improve the lepton identification efficiency. Such contributions are not included in 
this estimate. 



VII DO SEARCH FOR SUSY PAIR PRODUCTION FOLLOWED BY LQD^ 

LSP DECAYS 

S. Banerjee, N.K. Mondal, V.S. Narasimham, N. Parua 

A Physics motivation 

Recent interest in R-parity violating (RPV) SUSY decay modes is motivated by the possible high-Q^ event 
excess at HERA [33] . When interpretation of the excess through first-generation leptoquarks was excluded by 
the D0 [36] and CDF [37] experiments, it was suggested [34] that such an effect could be explained via the 
s-channel production of a charm or top squark decaying into the e+jet final state. Both the production and the 
decay vertices would thereby violate R-parity. Although more recent data has not confirmed the previous event 
excess, interest in RPV signatures has not abated. The CDF and D0 Collaborations have recently performed 
searches for RPV SUSY [25,26], and have set new mass limits on the RPV SUSY particles. Both experiments 
focussed their searches on the A' couplings, as motivated by the high-Q^ HERA event excess. The results of 
the D0 searches are here extended to the Run II case and the expected sensitivity to the RPV couplings is 
discussed. 
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FIGURE 23. Points in the {'mo,mi/2) SUGRA parameter space used to generate RPV events in the ee + 4 jets 
channel. 



B D0 Search for RPV neutrahno decays 



The D0 search for RPV SUSY considered the case of neutralino LSP which decays into a lepton and two 
quarks due to a finite RPV A' coupling (see Fig. 2). Both the electron and muon decay channels were considered, 
corresponding to what is commonly referred to as X'^j and X^ij couplings, respectively. The corresponding 
final states contain either 2e or 2^ and at least four accompanying jets. Unlike at HERA, this search is not 
sensitive to the value of the RPV coupling, as long as it is large enough so that the neutralino decays within 
the D0 detector. That corresponds to A' > 10^^, which gives a lot of room, given current indirect constraints 
(see Sect. HA). 



TABLE 5. Efficiency x BR (%), signal cross section and the expected 
event yield in 2 fb^^ of data, at various (ttio, mi/2) parameter space points. 
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TABLE 6. Expected lower limits on the squark and gluino masses for the case of no SUSY 
detection at Run 2. 
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We assume that the neutralino (LSP) pairs arc produced in cascade decays of other supersynimetric particles 
and use all SUSY pair production mechanisms when generating signal events. 

Signal events were generated within the SUGRA framework with the following values of SUSY parameters: 
Ao — 0, n < and tan/3 = 2 (the results are not sensitive to the value of Aq .) Center of mass energy of 
the colliding beams was taken to be 2 TeV. ISAJET [53] was used for event generation. The acceptance and 
resolution of the D0 detector were parameterized using the following resolutions: SE/E = 2%®15%/v'^ [GeV] 
(electrons), 5{l/p)/il/p) = 0.018 0.008(l/p) (muons), and SE/E = 3% ® 80%/Ve [GeV] (jets) and found 
consistent with the full detector simulation based on GEANT [69] . 

Figure 23 shows the points in the {rno, 1711^2) SUGRA parameter space where signal Monte Carlo events were 
generated for the electron channel. Similar points were studied for the muon-decay channel. 



C Selection criteria for the dielectron channel 

A multijet trigger was used for the analysis of Run 1 data. It was found to be nearly 100% efficient for the 
typical RPV signal. Since the Run 2 trigger list will include a similar trigger, we assume trigger efficiency of 
100% and do not perform any trigger simulations for the Run 2 analysis. 

The following offline selections were used: 

• At least two good electrons, the leading one with Erie) > 15 GeV and the other one with Erie) > 10 
GeV; 

• Rapidity range | ?7 |< 1.1 (central calorimeter), or 1.5 <| r/ |< 2.5 (end calorimeters) for all the electrons; 

• Energy isolation for the electrons: the EM energy in the R=0.2 cone about the center of gravity of the 
EM cluster, subtracted from the total energy in R=0.4 cone, should not exceed 15% of the EM energy 
in the R-0.2 cone. 



At least four jets with Et{j) > 15 GeV and \ rj \< 2.5; 
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FIGURE 24. Estimated exclusion contour for Run 2 in the (mo,mi/2) plane for tanl3 = 2, Aq — 0, /j. < 0, from the 
ee + 4 jets channel, assuming no SUSY signal is observed. Scenario I corresponds to a background of 36 ± 4 ± 6 events 
(direct scaling from Run 1); scenario II uses the background of 15 ± 1.5 ± 1.5 events (scaling, but with improvements 
in the detector taken into account). 



• The dielectron invariant mass (Mge) should not be in the Z-mass interval, ie, | M^e — Mz \> 15 GeV/c^. 

In the present analysis we have dropped the requirement on Ht = "^Erie) + "^Exij) , but retained all 
other offline criteria that were used in the previous analysis of data from Run I [26] . 



D Selection in the dimuon channel 

The following event selection requirements were used for the muon decay channel: 

• Two muons, the leading one with px > 15 GeV, and the other one with pT > 10 GeV. 

• Rapidity range |?7| < 2.3 for both muons. 

• Energy isolation requirement for both muons, i.e. the calorimeter energy accompanying the muon in a 
(rj (j)) cone of 0.4 should be consistent with that from a minimum ionising particle. 



At least four jets with Et{j) > 15 GeV and \r]\ < 2.5; 



R — parity violating SUSY in the jets + muons channe 




FIGURE 25. Estimated exclusion contour for Run 2 in the {mo,mi/2) plane for tanl3 = 2, Aq — 0, /j. < 0, from the 
fifi + 4 jets channel for background of 10 ± 1.0 ± 1.0 (direct scaling from Run 1). 



E Signal efRciencies 



The number of signal events expected can be written as: (N) — C ■ a ■ e, where (N) is the expected number 
of events for luminosity C, a is the cross-section, and e is the overall efficiency. The efficiency e can be split 
into three terms: e = etrig • fkin • Cid- Here etrig is the trigger efficiency for the events that pass the offline cuts ( 
assumed to be 100%), Ckin is the efficiency for offiine criteria, which includes kinematic, fiducial and topological 
requirements, and ej^ is the electron/jet identification efficiency. 

The efficiency for identifying jets is very high (> 95%) and is expected to stay the same in Run 2. 

Electron identification efficiencies in Run 1 were 80±7% in the central (\ri\ < 1.1) and 71±7% in the forward 
(1.5 < \ri\ < 2.5) regions [26]. These efficiencies were calculated for electrons with Exie) > 25 GcV, ft drops 
by about 30% for electrons with Erie) ^ 10 GeV. 

The muon identification efficiencies used in Run 1 were 62 ±2% in the central (|?7| < 1.0) and 24 ±4% in the 
forward (1.0 < jryj < 1.7) regions [70]. These were calculated for muons with pT > f5 GeV. For muons with 
10 GeV < pt < f 5 GeV the efficiencies were 80% smaller on average [71]. 

In the present analysis we have taken the overall particle identification efficiency to be 0.90 ± 0.09 in each 
channel, independent of lepton Et, primarily due to the expectation of a better tracker and muon spectrometer 
for the upgraded D0 experiment. 



F Backgrounds 

The main backgrounds are expected to arise from DrcU-Yan production in association with four or more jets, 
dilepton top-quark events, and QCD multijet events. The latter is the dominant background for the electron 
channel (followed by the Drell-Yan background). In the case of muons, the background is dominated by the 
Drell-Yan and top pair production. We used Monte Carlo to calculate background from the first two sources, 
and data to estimate background from QCD jets. 

Background for the Run 1 analysis was estimated to be 1.8 ± 0.2 ± 0.3 (with 1.27 ± 0.24 from QCD and 
0.42 ± 0.15 ± 0.16 from the other processes) for ^ 100 pb~^ of data. To extrapolate this number to the data 
set from Run 2, we have simply multiplied it by the ratio of luminosities to obtain 36 ± 4 ± 6 events. However, 
it is expected that due to the central magnetic field in the upgraded D0 detector, the probability of jets to be 
misidentified as electrons will be reduced by a factor of ~ 2 in Run 2. We have therefore considered a second 
scenario with the smaller expected background of 15 ± 1.5 ± 1.5 events. 

For the muon channel, the expected background has been scaled directly from the Run 1 analysis. We expect 
10 ± 1 ± 1 background events in Run 2. 

G Results 

In order to obtain the sensitivity of Run 2 in to RPV decays, we calculated the efficiency for signal for all the 
mass points shown in Fig. 23. Typical efficiencies, the signal cross section in the ee + 4 jets channel, and the 
expected event yield in 2 fb^^ of data, for several representative (mo, mi/2) points, are given in Table VII A. 
Similar numbers are obtained for the muon channel. 

We use these efficiencies to obtain exclusion limits in the (^0,7711/2) plane at 95% CL, assuming that no 
excess of events will be observed above the predicted background. The exclusion contours for the electron and 
muon channel are shown in Fig. 24 and 25, respectively. Numerical values of the limits are summarized in 
Table VII A. 

It's worth mentioning that our analysis provides a conservative estimate of the sensitivity achievable in 
Run 2, since no formal optimization of the signal vs. background has been performed. We expect that a formal 
optimization can improve the sensitivity in the mass reach by 15-20%. 

VIII SPARTICLE PAIR PRODUCTION WITH BARYONIC DECAYS: 
DIFFICULT SCENARIOS FOR SUSY DETECTION AT THE TEVATRON 

H. Baer, K. Cheung, J. Gunion 

A Motivation 

GUT scale boundary conditions {e.g. the niSUGRA universal boundary conditions) in which the gaugino 
soft-SUSY-breaking masses are assumed to have a common value at Mu lead to ratios of 7:2:1 (at low energy) 
for m-g : 'm~± : m^o due to M3 : M2 '■ Mi = a^ : a2 '■ cxi. In particular, the LSP Xi is bino-like and thus 

substantially lighter than the wino-like Xi- However, there are well motivated models in which Mi > M2. 

• In the context of string models, if the dilaton F component does not acquire a large vacuum expectation 
value (the most natural result in many models), then the leading contributions to gaugino masses will 
arise from loop corrections. This was first noted in the context of the the O-II model of Ref. [72] in the 
limit where supersymmetry breaking is ascribed entirely to the overall size modulus with none arising 
from the dilaton field. Including the possibility of a Green-Schwarz mixing term, parameterized by SqSj 
one has (at Mu): 

M3 : M2 : Ml °^'' -(3 + 5gs) ■ (1 - Sas) ■ (^ - Sas) ■ (18) 

5 

At low energies, one finds M2 : Mi = 2(1 — 6gs) ■ (33/5 — Sas)- For das in the preferred range from to 
—4, M2 < Ml and both the Xi and Xi are wino-like and very degenerate. More recently [73,74] the above 



boundary conditions have been re-discovered (at least in the Sgs = case) in a more general context. 
If the soft masses are dominated by the conformal anomaly, one obtains (to all orders in perturbation 

theory) Mi = '^ a m3/2- Going to low energy scales this translates into 

11a a 3as ,_,^, 

^1 = 1 2a~^3/2 , M2 = . 2^ "^3/2 , Ms = ^—^3/2 , (19) 

as consistent with the evolved result of Eq. (18) when Sqs ~ 0. 

• M2 < Ml at low energies can also arise if SUSY is broken by an F-term that is not an SU(5) singlet [75]. 
Recall that the Ma arise from a chiral superfield <& that appears linearly in the gauge kinetic function 
with (F) 7^ for its auxiliary component. 

jr = J d^eWaW,^ + h.c. ~ ^^A,A, , (20) 

where Xa,b {a,b = 1,2,3) arc the gaugino fields. In the general case, $,i^$ G (24x24)symmotric = 
1 ® 24 75 © 200. Only F$ components that are 'neutral' under SU(3), SU(2), U(l) can have a vev if 
these groups are unbroken after SUSY breaking. As a result, {F,p)ab — CaSabj with Ca depending on the 
representation. Ca is independent of a for the 1 representation (universal Ma at Mu) but M2 < Mi for 
the 200 representation. 

To give some specific resuhs, at Mu (Mq) the O-II model with Sgs = -4 yields M 3 : M 2 : Mi = 1 : 5 : 53/5 
(6 : 10 : 53/5) while the 200 model yields 1 : 2 : 10 (6 : 4 : 10). Even though M2 is only slightly less than Mi 
at low energies in the O-II and 200 models, Am<^ = 'm~± — ttijo < a few GeV is typical. Note that in both 

models nig is close to the common Xi 7 Xi mass. 

The phenomenology for a model in which M2 < Mi and R-parity is conserved was explored in [76,77] 
(see also [78]). There it was noted that at an e+e^ collider it will be virtually impossible to see Xi'xr P^i'" 
production since the potentially visible particles in Xj^ ^ Xi + • • • decays will be very soft. One can evade this 
problem by tagging the Xi'xr event with a photon, i.e. look for e+e^ — > 7Xi'xr- Sensitivity is substantially 
reduced compared to normal mSUGRA boundary conditions. At hadron colliders, if 'ni~± ~ m^o one loses the 

tri-lepton signal from X2X1 ~^ i^i~i'^ +i/t since the i'^ is too soft to be detected. Further, the like-sign 
di-lepton signal from gg production (for instance) deriving from gg -^ 4j -I- Xi Xi ~^ 4j -I- £ £ +i/t will also 
be undetectable. Still, as summarized in more detail later, in an R-parity conserving model the iets+ j^t signal 
will remain and allow reasonable (although somewhat reduced compared to mSUGRA) sensitivity to SUSY 
particle production [77]. 

B Baryon Number Violation 

This leads us to the worst-case scenario in which we allow for baryonic R-parity violation. We consider 
the case where the leptonic R-parity violating terms are zero (A^fe, A^-^., Ki — 0) and only A" is non-zero. We 
note that A''„„ must be antisymmetric in the n, m family indices (due to the implicit antisymmetric colour 
coupling). For example, converting to a notation where p, n, m are given in terms of the quark type, we might 
imagine A^;^^ = -A^;^^^ ^ 0. ^ 

The most obvious implication of baryonic R-parity violation is the possible absence of any specifically super- 
symmetric source of missing energy. The decay products of the Xi-LSP^ would be jets that would be visible in 
the detector provided the decays are prompt. Combining baryonic R-parity violation with Xi — Xi degeneracy 
provides some very unusual and, in some cases, very difficult scenarios for detecting supersymmetry. 

Of course, there are various means for placing bounds on some of the A"-^, (see, e.g. , [79] and [2]). Pertur- 
bativity up to the GUT scale requires that all A"'s be < 1. The experimental limit on the AB — A5 — —2 



decay ^^O -^^^C+K++K+ [14,80] leads to |A;:,,| < 5.6 x 10"^ ('-f^) [t^B^Y" {t4v) ■ Here, 



^' We do not have space here to consider the p-LSP case that is also possible for O-II model boundary conditions when 
5gs ~ —3. 



^QCD is a scale that dimensional analysis suggests should be of order Aqcd ~ 200 MeV. Constraints on 
^ubd' ^tsd ^^"^ Kbd from N — > N osciUation limits (see Ref. [79]) vary roughly as l/Mj^j^, where M^j^ 
is the L — R mixing entry in the relevant squark mass-squared matrix. We assume that L — R mixing can 
be neglected, implying that Af^_^ is small and that these bounds are not relevant. 2^2 processes such as 

u + d —> d* ^f dxl disturb the standard baryogenesis picture if sphalerons are present and if any |A"| were 
larger than ^ 5 x 10^^(?7i^/l TcV)^/^, a very severe constraint. However, these constraints do not apply if 
sphalerons are not present at the appropriate epoch [8] . Bounds that can be placed on various products of the 
A"'s are not relevant to our analysis since we shall consider only one coupling to be non-zero at a time. 



C Phenomenology 



X° decays 

The expected ct for the Xi is of considerable importance. If the Xi ~^ Ufjdnd'j^ (family indices implied) 
decay is prompt, discovering supersymmetry at a hadronic collider in this model is far from easy. However, if 
this decay has a significant path length, our ability to search for supersymmetry would be greatly enhanced by 
virtue of being able to look for a decay vertex. In the limit of a very long path length, we would revert to the 
jets+i/y signature. This was studied in Refs. [76,77] for the present class of models. Since m 



~xt 



m^a and my. 



is not that much larger than ra~± supersymmetry signals are weaker than for mSUGRA boundary conditions. 

Still, the LHC would be guaranteed to find a signal for m-g ^ 1 TeV; but at the Tevatron the weaker signals 
would imply a fairly limited discovery reach. 

A result for the path length is easily given in the limit where the sfermions are much heavier than the 
x\- Assuming also that the relevant w-type and d-type squarks are more or less degenerate with mass m and 
neglecting L — R squark mixing, one obtains 
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(21) 



where Nij are the neutralino mixing matrix entries with i refering to Xi'=i 234 ^^.^ J refering to -B, 14^3, -ffj*, H2- 
Eq. (21) neglects the tiny contributions from the Higgsino components of x"- Iii the absence oi L — R squark 
mixing there is no contribution to the RPV decay of the Xi coming from its primary wino component (of size 
A^i2). (This is because the RPV A" superpotential term only involves right-handed superfields.) In the model 
class being considered, Nn ^ 0.1 is typical and fh > 1 TeV is preferred (although this could be avoided in a 
'no-scale' type model). For to -- 1 TeV, A" < 0.017 (< 0.0043) is required for ct > 100 iim if m^o ~ 100 GeF 
(^ TO-i), respectively. For m ^ 300 GeT^, the corresponding A" values must be smaller by a factor of .09 for 
CT > 100 iim. Given that typical Yukawa couplings for the heavier 'light' quarks are gms/{2mw) ^ 0.0008 
and gmc/{2mw) ^ 0.006, the RPV analogues could easily be this small. 

For the constraint given earlier on A"^^, the Xi would have a very long path length if all other A"'s are zero. 
However, it seems likely that the higher generation A"'s would be larger than A"^^. The largest A"'s might 
be those associated with the 3rd family. However, since the mass of the Xi is typically not very large in the 
models we consider, a top quark in the final state might have to be virtual, implying substantial phase-space 
suppression. The values for ct x [|A"|^|A^ii|^] for t'^*^bd^ t'^*^hs and t^*^sd are plotted in Fig. 26 in comparison 
to that obtained in any one of the usd, ubd, ubs, csd, chd and chs channels. At m^a — 100 GeV^ (175 GeV^), 
one of the three virtual top channels would dominate over one of the other six channels if A" for the former 
is a factor of ~ 5 x 10** (^ 500) larger than for the latter. Keeping in mind that Yukawa couplings differ by 
ratios of ^ mt/niu ^ 5 x 10^ and ~ mt/rng ~ 350, such large ratios should not be regarded as impossible. 

Typical ct values for the various decays are obtained from Fig. 26 by multiplying by l/|A'^ii|^ ~ 100 and by 
1/1 A" 1^. If Xi decay is dominated by a virtual top channel and the relevant A" ~ 1, a path length of order 
100 m (100 fim) is to be expected for m^o ^ 100 GeF (~ mj). For smaller A", ct is even larger. However, if 



one of the Xusd, 
xt decays 



is substantial, the Xi decay will be very prompt. 



In the models being considered, Am^ — m~± 



< 1 GeV is typical and the ct for Xi 



X?. 



(see 



[78] ) can be so large that the ct for direct RPV decay of the Xi (values being comparable to those plotted for 
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FIGURE 26. We plot ct x [|A"| lA'^ii] ] as a function of m^o for a variety of final states, including those arising from 
virtual top intermediate states. 



the Xi) can easily be shorter (i.e. the width larger). In such scenarios, direct Xi decays to the xl rnight be 
quite rare. Phenomenologically, however, it would hardly make a dilFerence since the decay products of the Xi 
decay are so soft for the small expected values of Am^ anyway. 
g decays 

The normal SUSY g -^ qqxl amd g -^ q'qXi decays compete with the RPV g — > qq'^ — ^ 3q decays. For 
our boundary conditions, the former SUSY modes can be sufficiently phase space suppressed that the RPV 
decays could dominate if one or more of A"^^, A"{,^, A"^^, A"j,^, A"^^, or A"^^ were large. However, this would 
result in little change in the phenomenology since the extra 2 jets coming from the SUSY g chain decays would 
generally be quite soft and not easily reconstructed at a hadron collider. At an e~^e~ collider, g production 
(via e+e" -^ qqgg) has a small cross section and g decays are not an important issue. 



D SUSY detection at a hadron collider 



Various possible scenarios could emerge. Suppose first that the relevant A" is sufficiently small that the Xi 
is effectively stable in the detector. Then, the R-parity conserving results of Ref. [77] for nig ^ "^y± — "^i^o 

apply. Since the lepton from Xi ~^ ^ Xi decay is very soft, standard SUSY signals relying on leptons are no 
longer viable. This includes the like-sign dilepton signal for gg production in which gg — > 4j + Xi Xi + X -^ 
4j + ^+£+ +-^T (or the charge conjugate) and the trilepton signal from Xi X2 ~^ l^^i'^i'" . One must rely 
entirely on the jets plus missing energy signal for gg production, which is weakened by the fact that the jets in 
both g — > q'qxi ^^d Xi ~> q'qXi decays are no longer hard enough to pass typical cuts used by CDF and DO to 
define a visible jet. The only energetic high-p^ jets present are those arising from initial state radiation. Due 
to the smaller number of jets passing cuts, discovery reach at the Tevatron is substantially reduced compared 
to the normal mSUGRA scenario, even if the relatively soft (as appropriate for the scenario being considered) 
Run I CDF and DO cuts are employed. This is illustrated in Table 7. 

Next, suppose A" is of 'moderate' magnitude such that Xi decay is visible in the detector, i.e. has a path 
length leading to a visible vertex in the inner vertex detector or tracker or a recognizable delayed decay in the 
outer portions of the detector. The combined requirement of many jets and a clear vertex or delayed decay 
should leave very little QCD background. Supersymmetry discovery would be largely limited by event rate 
and discovery reach would exceed even that typical of mSUGRA boundary conditions. 

The most difficult case for supersymmetry discovery would be that in which A" is sufficently large that the 
RPV Xl decay is prompt and the A"'s involving the top/stop are sufficiently small that the virtual or real 
top-quark decay modes of the Xi are not important and we must consider x^ — > 3j decays. In this case, 



TABLE 7. Maximum mg values that can be probed using DO and CDF Run I cuts 
in the jets +IfT final state for different integrated luminosities, L, at the Tevatron 
and Tev33 using the O-II model boundary conditions with Sqs = —4.5. Observ- 
ability is defined by S/y/B > 5 and S/B > 0.2. Also given are the maximum mg 
values for which gg production can be observed in the case of mSUGRA universal 
boundary conditions for both Run I cuts and also stronger cuts optimized for Run 
II/TeV33. The results of this table are for tan/3 — 2 and fi < 0. 





Cuts 


L = 


8 pb^' 19 pb-^ 100 pb-^ 2 fb"^ 25 fb"^ 


DO 
RunI 


Sgs = —4.5 
mSUGRA 


— SOGeV 110 Gel/ 150 GeV 150 GeV 

— UOGeV 200 GeV 200 GeV 200 GeV 


CDF 
RunI 


Sgs = -4.5 
mSUGRA 


100 GeV 130 GeV UOGeV 160 GeV" 160 GeV 
150 GeV 170 GeV 210 GeV 240 GeV 240 GeV 


Strong 


mSUGRA 


rUg excluded by m^± < 47 GeV 250 GeV 300 GeV 



3j = uds, udh, usb, cds, cdb, and/or csb are all, in principle, possible. There would be no specifically SUSY 
source of missing energy in sparticle-pair production. The most characteristic feature of gg production would 
be the presence of many energetic jets in most events. To explore the characteristics of such events in more 
detail, we have considered a sample case defined by the following O-II (GUT scale) model parameters: 



leading to 



^GS 



A2 



-4.1, 



mo = 1800 GeV, mi/2 == 10 GeV, Aq ^ 0, tan /3 = 3, M > , 



102 GeV, 



m^± 



100 GeV, At 



0.23 GeV, 



771 c;o 

X2 



120 GeV, 



-^3 



771^0 
A4 



1140 GeV, 






1.4- 1.8 TeV. 



(22) 



(23) 



Note the very low mass scales for all gauginos, including the g, and the high mass scales (typical for the O-II 
boundary conditions) of all sfcrmions and higgsinos. 

We have performed a Monte Carlo simulation of signal and background (using IS A JET [53]). We define 
jets using Et > 15 GeV, \r]\ < 4, and a cone size of AR = 0.5. In Fig. 27(a), we plot the event rate for the 
QCD background as compared to that for supersymmetric particle pair production (all channels included, but 
mainly gg) as a function of the number of jets, rij, before any tagging requirement, and before any other cuts. 
The signal has larger Uj on average, but the QCD background is still bigger by a factor of 10 even if we require 
Uj > 8. Nonetheless, if we compute the nominal S/^/B for the rij — 8 bin, assuming L = 2 fb^ (30 fb^ ), 



we get S/VB = 45 (170). Thus, if the systematic uncertainty for the rij = 8 background could be reduced 
to below 10% we would achieve an observable signal without tagging. This might be possible by normalizing 
using the rij < 7 bins once a reliable calculation of relative rates for multi-jet events with rij ranging from low 
values up to rij = 8 becomes available. (The ISAJET shower model results for relative rates are probably not 
accurate at the 10% level.) 

However, the very strong limit on A"^,^^ means that one (or more) of the ubd, ubs, csd, cbd, cbs channels will 
dominate x? decay. All of these contain at least one heavy quark. Since SUSY pair production events contain 
at least two Xi's, we can then use double-tagging to reduce the background relative to the signal. For all but 
the cds decay channel, at least two of the jets in the final state will be 6-jets. If cds decays dominate, we would 
have two c-jets. Since c-tagging will have lower efhciency than 6-tagging, dominance of xl decays by the cds 
channel results in the most difficult scenario for supcrsymmetry detection. 

To explore this 'worst case', we assume Xi ~^ csd exclusively. Fig. 27(b) gives the rij distributions for 
signal and background after requiring that every accepted event contain at least two tagged jets. For this 
plot we use the CDF 'JP'-tag algorithm (designed for c quark tagging). We require that a tagged jet have 
|r/| < 2 '' and pr > 25 GeV. For such a jet, the JP-tag yields [81]: et > 0.5, Ec > 0.3, for €q,g - 0.05, 



where q 



u, 



From Fig. 27(b), we see that after tagging it is best to focus on the rij — 7 bin, with 



*' This is the rapidity coverage for Run II. 
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FIGURE 27. The number of events is plotted as a function of the number of jets for the QCD background (dashes) 
and for SUSY production processes (sohd) (a) without tagging (b) after requiring that at least two jets be tagged (using 
the CDF JP charm tagging algorithm). O-II model parameters as in Eq. (22). 



S =lAx 10'', B == 3.2 X 10'' (5* = 2.1 x 10^ B = A.8x 10^) for i = 2 fb"^ (L = 30 fb"^), leading to S/y/B = 78 
{S/vB — 303) with S/B ~ 45% (a big improvement). Thus, for low- mass SUSY scenarios [see Eq. (23)] of 
the type being considered, double-tagging leads to a very clear excess (beyond Standard Model expectations) 
of events with a large number of jets. 

At the low rrig value considered above, additional cuts in the variables E^* = ^ ■ \Ej,\ and Ht = X)? I-^tI+^t 
are not particularly helpful. However, in order to reduce the QCD background when searching for gluinos of 
substantially higher mass, cuts requiring high values for one or both of these variables are essential. For an 
appropriate cut, we would again examine the rij distribution after requiring two tagged jets. We anticipate 
that a signal will be observable so long as the basic production rate for gg pairs is substantial. The exact mass 
reach (in rrig) is under study [82]. 



E Discussion and Conclusions 



We have shown that observation of a signal from gg pair production is likely to be possible at the Tevatron 
even in the extreme of O-II model boundary conditions and RPV 3j decays of the Xi^ so long as the gluino mass 
is not too large. However, the high-jct-number signal upon which we focused does not allow us to determine 
the strength of the dominant RPV A" coupling. There are two possibilities for directly determining A". If A" 
is large, RPV-induced single squark production cross sections are also typically substantial, and, if a signal 
can be isolated [82], the cross section size gives a measure of A". If A" is small, the decay path length for the 
Xi might be observable and would again provide a measure of A". However, our preliminary estimates [82] 
indicate that there is a region of intermediate A" at higher 7715^ for which neither ct(xi) nor a(qfi) will be 
measurable. In this region, determination of A" would only be possible if an RPV decay mode of the Xi is 
competitive with its standard SUSY decay modes and these can be separated from one another (which is not 
at all certain); the relative size of the branching ratios would then provide a measure of A". Parenthetically, 
we note that at an e+e~ collider the only means for measuring A" would be via the Xi decay length; there are 
no sources of quarks or antiquarks as needed for squark production via baryonic RPV couplings. 

To summarize, even if there is baryonic R-parity violation, implying absence of missing energy in SUSY 
events, the very large number of jets expected for a typical SUSY event will allow us to isolate the SUSY signal 



so long as gluino (and/or squark) masses are sufficiently low that SUSY event rates are substantial. Further, 
there is at least a decent chance that we will be able to measure the strength of the RPV coupling. 
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PART 3: Spontaneous R-Parity Breaking 

IX DECAYS OF THE TOP-QUARK AND THE TOP-SQUARK 

F. de Campos, M.A. Diaz, 0. J. P. EboU, M. B. Magro, L. Navarro, W. Porod, D. A. Re- 
strepo, and J. W.F. Valle 

A Introduction 

In this section we study the case of spontaneous R-parity breaking. This is a modification of the previous studies 
in that an extra superfields is added: vr, which is a gauge singlet and corresponds to a right-handed neutrino. 
Here we focus on the case of bilinear R-Parity Violation (BRPV) [83-86], where A, A', A" = in Eq.(3) and 
Ki ^ 0^. This includes as an additional feature a vev for the sneutrinos. These models are well- motivated 
theoretically as they arise as effective truncations of models where R-Parity is broken spontaneously [87] 
through right handed sneutrino vacuum expectation values (vev) {v^) ~ vj^^ 0. They open new possibilities 
for the study of the unification of the Yukawa couplings [88] . In particular it has been shown that in BRPV 
models bottom-tau unification may be achieved at any value of tan/3. From a phenomenological point of view 
these models predict a plethora of novel processes [89] that could reveal the existence of SUSY in a totally 
different way, not only through the usual missing momentum signature as predicted by the MSSM. They 
provide a very predictive approach to the violation of R-Parity, which renders the systematic study of R-parity 
violating physics [89] possible. 

We consider the simplest superpotential which violates R-Parity 

W"|{ = Wmssm + ^iL^Hu , (24) 

assuming that TRPV terms are absent or suppressed, as would be the case if their origin is gravitational [90]. 
The €i terms violate lepton number in the ith generation, respectively. Models where R-Parity is broken 
spontaneously [87] through a vev of the right handed sneutrinos (i>^) = vr ^ Q generate only BRPV terms. 
The Ci parameters are then identified as a product of a Yukawa coupling and vr. This provides the main 
theoretical motivation for introducing explicitly BRPV in the MSSM superpotential. For simplicity we set 
from now on ei = €2 = 0. Thus we have tau-lepton signatures and we only test the case where tau-lepton 
number is violated. The MSSM-BRPV has the following superpotential 



Wfi,^ Sab 



htQlU^Hl + hbQlD:iH2 + hrL\R^H2 + iiKH''^ + e:iLlHl , (25) 



where the first four terms correspond to the MSSM. The last term violates tau-lepton number as well as 
R-Parity. 

In this model the top-quark and the top-squark get additional decay modes, e.g. t ^ fj^ 6 or ti —tr^+h. We 
study these decays in view of the Tevatron (for top decays in TRPV models see [91]). We show that existing 
Tevatron data give additional constraints on the parameter space and extrapolate to Run II. 

B Top Decays 

One of the major successes of the Tevatron has been the discovery of the top-quark [92]. The large top mass 
implies a relatively small production cross section at the Tevatron. Therefore, the sum of all branching ratios 
of the top decays except t — > W^ b is only restricted to be smaller than approximately 25 % [93]. In the MSSM 



Throughout the rest of this section we shall denot the parameter k as e. 
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FIGURE 28. (a) Branching Ratios for f — ^ 61 r"*" as a function of £3. (b) Branching Ratios for t ^ Tj^ fe as a function 
of W3 for different values of tan/3. In both cases the parameters are: M = 180 GeV, fi = 200 GeV, tan/3 — 35, 
Me^ = 285 GeV, A^ = 280 GeV, Mq = 285 GeV, Mu = 180 GeV, Md = 190 GeV, At = 320 GeV, At = 120 GeV, 
B = 50 GeV, -500 GeV < e < GeV, GeV < B2 < 200 GeV, 1 GeV < 1^3 < 50 GeV. 

the top can decay according to: t -^ W~^ b, t ^ H^ b, t ^ Xiii, t ^ Xi ■ "^ In the BRPV model the charginos 
mix with the charged leptons, the neutrahnos with neutrinos, and the charged sleptons with the charged Higgs 
boson [83-85]. Therefore, the top can have additional decay modes: 



t 



t —t VtIi 



t 



bi. 



(26) 



As an illustrative example, we show in Fig. 28a the branching ratio for i — > t+ &i as a function of a) 63 and b) 
W3. We have randomly chosen 10,000 points imposing the following experimental constraints: m^^ < 18 MeV, 



"ii' 



> 80 GeV, ■min{mu+,mf^) > 70 GeV, and 



"^ict 



> 85 GeV. The parameters are listed in the 



figure caption. We find a strong correlation between the R-parity decay branching ratio BR{t 



h) 



and the magnitude of £3. This can be understood in the following way: in the chargino mass matrix the 
mixing between the leptons and the charginos disappears if one does the following rotation in the superfields: 
Hi -^ N{^Hi — £3^3) and L3 -^ N{ijlL^ + e^Hi) (N being the normalization). In this basis the coupling 
between i, r, and 61 is proportional N hb £3 leading to this feature. 

In Fig. 28b wc show the branching ratio for t ^ f^ b a.s a. function of W3. As in Fig. 28a, results in Fig. 28b 
are displayed for different values of tan /3 and the other parameters are also the same. The dependence on tan/? 
is a result of: (i) The stau - charged Higgs boson mixing is proportional to the R-parity breaking parameters 
£3 and U3 (ii) The decay width depends on the bottom Yukawa coupling which increases with tan /5. As can be 
seen from the figure, there is a strong correlation between the magnitude of the R-parity breaking branching 
ratios and the mixing between the stau and the charged Higgs boson. 

We have performed a similar scan for small tan/3 for both of BRPV decay channels discussed above. These 
are suppressed in this case and can not exceed 2% or so, i.e. ^ BR{t -^ bX) < 1 — 2% {X ^ W)), because their 
decay widths are in all cases proportional to the bottom Yukawa coupling squared. In the case of t — > r^ bi 
this is clear from the discussion of Fig. 28a. In the case of i ^- f^*" 6 one has to note that the stau mixes mainly 
with the charged component of the down-type Higgs multiplet Hi {jl and Hi have the same gauge quantum 
numbers) and the Hitb coupling is proportional to hi,. 

In every case the various decay modes lead to the cascade decays 
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(27) 



The mode 61 is practically ruled out by existing LEP2-data [94]. 



In nearly all cases there are two r's and two 6-quarks in the final state plus the possibility of additional leptons 
and/or jets. Therefore, 6-tagging and a good t identification are important for extracting these final states. 
Moreover there is in general a large multiplicity of charged particles in the final state which should be helpful 
in reducing the background. The background will come mainly from the production of one or two gauge bosons 
plus additional jets. The conclusion in similar cases [95] has been that in its next run the Tevatron should 
be sensitive to branching ratio values up to 10"'^ — 10^^ depending on the mode. Therefore, the possible 
observation of one of these additional decay modes at the Run II of the Tevatron should give a strong hint on 
the underlying parameters. 



C Top-Squark Decays 

Top-squark physics is a very interesting part of supersymmetric theories, because the lighter top-squark might 
be the lightest charged SUSY particle [58]. In the kinematic region accessible to the Tevatron the light top- 
squark has the following MSSM decay modes: h ^> xt ~^ ^> ^i ^^ Xi + c, h -^ if + vi + b, ii ^ vi + l^ + b, 
^1 ^^ Xi + W^~^ + ^: ^^^d ii ~^ Xi + H^ + b (for a discussion see e.g. [96] and references therein). In BRPV 
models the top-squark has an additional and phcnomcnologically very interesting decay mode [97]: 



ti^T+ + b (28) 



In the following we have concentrated on scenarios where only the two-body decay modes are possible. We adopt 
the framework of Supergravity unification as in [83] in order to reduce the number of free SUSY parameters. 
However, we keep es and V3 as free parameters for the moment. In Fig. 29a we show the areas in the TOj^-to^o 
plane where the branching ratio ti ^ t+ + 6 is larger than 90% for different values of £3 and v^. We restrict 
to the range [eal, {v^l < 1 GeV, and vary randomly the MSSM parameters keeping rrif < m~± + nib. This 
demonstrates that one can get a dominance of the R-Parity violating decay mode even for relatively small 
values of the R-parity breaking parameters. The upper-left triangular region corresponds to ttij: < m^a + jn^ 

and thus BR{ti -^ br) = 1. In the lower-right triangular region m^ > m~+ + nii, and therefore ii -^ bxt is 

open. In the central region the top-squark has the two decay modes ii —> br and ii —> cxi- The solid lines, 
defined by the maximum value of [esl and |w3|, are the boundary of the regions where BRiti -^ br) > 0.9 such 
that points at the left of the boundary satisfy that condition. 

Since BRPV models allow the decay {ti -^ t^ +b) we can interpret the top squark as a third generation 
leptoquark. Therefore we can use the limits obtained from leptoquark searches [36,37] to derive limits on the 
top-squark for this case. In Fig. 29b we show an exclusion plot in the mQ-mii2 plane. The nearly horizontal 
dashed lines are chargino mass contours and the lines forming radial patterns are the top-squark mass contours. 
The upper to the lower radial curves corresponds to mi_^ = 120, 100 and 80. The region limited by the dotted- 
dashed line is defined by mi_^ < rn-+ . The analysis rules out too and TO.1/2 points in the dark hashed region. In 
the lower hashed region no points with radiative electroweak symmetry breaking can be found. We have taken 
tan/3 ~ 3, Aq = —650 GeV and e^/fi = —0.5 and verified that in this region BR{ii -^ br) = 1. The Tevatron 
limits can not be directly applied when Aq > —500 GeV, because in this case rrif > m~+ + to,;,. The regions 

1 Xi 

in the too-TOi/2 plane where m^ < to,^+ are excluded if —650 < Aq < —500 GeV and |e3//i| is sufficiently 
large so that the three-body decays are negligible. Fig. 29a shows that the region where mz < m~+ + to.^ and 

1 Xi 

BRiti -^ T+ + fe) « 1 is practically ruled out by experiment. For this particular choice of SUSY parameters 
there is only a little window still to explore at the Run II of the Tevatron. However for other choices of SUSY 
parameters, e.g. Aq = —900 GeV the dark-hatched region fills up only about half of the allowed region where 
TO.J < m~+ + rrib and would therefore be open for investigation at the next run. 

Xi 

The MSSM three-body channels could be competitive with the BRPV one if |e3//i| is very small and mo ^ 
TO.i/2- In this case the condition BR{ti ^ t+ + fe) « 1 no longer holds and our analysis is not applicable. If 
ks/z^l < 10^^ (which leads to tau neutrino mass in the 10~^eV range in the mSUGRA model) for the same 
value of tan/3 we find [98] that the decay mode into cyfl is competitive with the br channel. The ii -^ cyfl 
channel becomes more important for large tan/3 and to^ < to~+. In this case one needs |e3//x| > 10^^ in 

order to get a neghgible BR{ii -^ ex"). 
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FIGURE 29. (a) Contour- lines for BR(fi -^ far) > 0.9 in the mj^-m^o plane. The gray region shows the area where 
only those two decay modes are open. We consider |e3|, \vz\ < 1 GeV, and the MSSM parameters are varied randomly 
such that mj^ < m.± + mt- The lines are defined by the maximum value of |e3| and \v3\ and delimit the regions where 

BR{ti -^ br) > 0.9. (b) Exclusion contour in the mo-mi/2 plane. The nearly horizontal dashed lines are chargino 
mass contours while the radial-like dashed lines are the top-squark mass contours. These change from solid to dashed 
when the top-squark becomes heavier than the lightest chargino. The radial curves correspond to m^^ = 120, 100 and 
80, respectively, from upper to the lower. The region limited by the dotted-dashed line has m^^ < ?n^+ . The dark 
hashed region is excluded by experimental data while the lower light-hashed region is disfavoured by theory. We have 
fixed tan/3 = 3, Ao = -650 GeV and es/m = -0-5. 



D Summary 

We have studied top-quark and top-squark decays in a supersymmetric model with bilinear R-parity break- 
ing. We have found that in both cases there exist additional top and stop decay modes leading to novel 
phenomenological implications with respect to those of the MSSM. In the top-quark case the new decay modes 
are t ^> fi h, t ^r Vrti, and t — > r^hi. We have shown that existing data on non-W top decay from the 
Tevatron are already sensitive to the BRPV parameters, adding both sbottom and stau decay channels. 

In this model the top-squark has the additional channel ii — > t+ + b. This channel will be 100% if the 
stop is the lightest SUSY particle, which is possible in the BRPV model. Moreover, we have demonstrated 
that this decay can be dominant even when the lightest neutralino below the stop and the R-parity breaking 
parameters \e^\ and \v^\ are well below a GeV, as long as the R-parity conserving chargino decay mode is 
kinematically closed, i.e. for mi < m -+ . We have studied scenarios in a SUGRA model with universality of 
the soft breaking terms at the unification scale and we have found that the Tevatron data on third generation 
lepto-quark data rule out scalar-top masses below 80-100 GeV, depending on the parameters. 
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Note Added in Proof: After this report was submitted a further slepton production process via R-parity 



violation was discussed in the literature [99]. 
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